Carleman estimates for semi-discrete parabolic 
operators with a discontinuous diffusion coefficient 
and application to controllability 

Thuy N.T. Nguyen* 
Abstract 

In the discrete setting of one-dimensional finite-differences we prove 
a Carleman estimate for a semi-discretization of the parabolic oper- 
ator dt — dx{cdx) where the diffusion coefficient c has a jump. As 
a consequence of this Carleman estimate, we deduce consistent null- 
controllability results for classes of semi-linear parabolic equations. 

1 Introduction and settings 

Let f],a; be connected non-empty open interval of M with a; (e il. We 
consider the following parabolic problem in (0, T) x Vl, with T > 0, 

dty - d^icdxv) = luiV in (0,r) x n, y\dn = 0, and y\t=o = yo, (1.1) 

where the diffusion coefficient c = c(x) > 0. 

System is said to be null controllable from yo £ in time T if 

there exists v € L'^{{0,T) x Q), such that y{T) = 0. 

In the continuous framework, we refer to |FI96) and |LR95j who proved 
such a controllability result by means of a global/local Carleman observabil- 
ity estimates in the case the diffusion coefficient c is smooth. The authors of 
|BDL07| produced this controllability result in the case of a discontinuous 
coefficient in the one-dimensional case later extended to arbitrary dimension 
by |LR10] . Additionally, a result of controllability in the case of a coefficient 
with bounded variation (BV) was shown in [FCZ021 IL07] . 

The authors of [LZ98j show that uniform controllability holds in the 
one-dimensional case with constant diffusion coefficient c and for a constant 
step size finite-difference scheme. Here, "uniform" is meant with respect 
to the discretization parameter h. The situation becomes more complex in 
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higher dimension. In fact, a counter-example to null-controllabihty due to 
O. Kavian is provided in |Zua06j for a finite-difference discretization scheme 
for the heat equation in a square. 

In recent works, by means of discrete Carleman estimate, the authors 
of [BHLlOa] . [BHLlOb] and |BL12| obtained weak observabihty inequahties 
in the case of a smooth diffusion coefficient c{x). Such observabihty esti- 
mates are charaterized by an additional term that vanishes exponentially 
fast. Morever, also with a constant diffusion coffiencient c, under the as- 
sumption that the discretized semigroup is uniformly analytic and that the 
degree of unboundedness of control operator is lower than 1/2, a uniform 
observability property of semi-discrete approximations for System (jl.ip is 
achieved in [LT06] . Besides that, such a result continues to hold even 
with the condition that the degree of unboundedness of control operator is 
greater than 1/2 [NT2] . 

In the case of a non-smooth coefficient, our aim is to investigate the uni- 
form controllability of System (jl.ip after discretization. It is well known that 
controllability and observability are dual aspects of the same problem. We 
shall therefore focus on uniform observability which is shown to hold when 
the observability constant of the finite dimensional approximation systems 
does not depend on the step-size h. 

In the present paper we prove a Carleman estimate for system (jl.ip in 
the case of: 

• the heat equation in one space dimension; 

• a piecewise coefficient c with jumps at a finite number of points in 
Q; 

• a finite-difference discretization in space. 

The main idea of the proof is combination of the derivation of a discrete 
Carleman estimate as in [BHLlOal IBL12j and tecniques of [BDLOTj for op- 
erators with discontinuous coefficients in the one-dimensional case. A similar 
question in n-dimensional case, n > 2, remains open, to our knowledge. 

When considering a discontinuous coefficient c the parabolic problem (jl.ip 
can be understood as a transmission problem. For instance, assume that c 
exhibits a jump at a E Jl. Then we write 

' dty - d^icd^y) = l^v in (0,r) x ((0, a) U (a, 1)) , 

' Cdxy\a+ = Cdxy\a-, y\a+ = Vla-, 

Mdn = 0, and y\t=o = yo- 

The second line is thus a transmission condition implying the continuity of 
the solution and of the fiux at x = a. 

When one gives a finite-difference version of this transmission problem, 
a similar condition can be given for the continuity of the solution. Yet, for 
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the flux, it is only achieved up to a consistent term. In what follows, in the 
finite-difference approximation, we shall in fact write 



(the discrete notation will be given below). Note that the flux condition 
converges to the continuous one if /i — > 0, /i being the discretization param- 
eter. This difference between the continuous and the discrete case will be 
the source of several technical points. 

An important point in the proof of Carleman estimate is the construction 
of a suitable weight function ip whose gradient does not vanish in the com- 
plement of the observation region. The weight function is chosen smooth in 
the case of a smooth diffusion coefficient c(x). In general, the technique to 
construct such a function is based on Morse functions (see some details in 
|FI96j ). In one space dimension, this construction is in fact straightforward. 
In the case of a discontinuous diffusion coefficient, authors of |BDL07j intro- 
duced an ad hoc transmission condition on the weight function: its derivative 
exhibits jumps at the singular points of the coefficient. In this paper, we 
construct a weight function based on these techniques in the one-dimentional 
discrete case. 

From the semi-discrete Carleman we obtain, we give an observability 
inequality for semi-discrete parabolic problems with potential. As compared 
to the result in continuous case |BDL07| the observability estimate we state 
here is weak because of an additional term that describes the obstruction to 
the null-controllability. This term is exponentially small in agreement with 
the results obtained in [BHLlOal IBHLlObj in the smooth coefficient case. A 
precise statement is given in Section [6l 

Finally, the observability inequality allows one to obtain controllability 
results for semi-discrete parabolic with semi-linear terms. In continuous 
case, this was achieved in |BDL07j . Taking advantage of one-dimensional 
situation, the results we state are uniform with respect to the discretization 
parameter h (see Section [6]). 

1.1 Discrete settings 

We restrict our analysis to one dimension in space. Let us consider 
the operator formally defined by ^ = —dx{cdx) on the open interval Q = 
(0,L) C M. We let a' G O and set Vti := (0,a') and := {a',L). The 
diffusion coefficient c is assumed to be piecewise regular such that 




n+V 



< c, 



< c < c, 
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with a € C'^{Vti),i = 1,2. 

The domain of A is D{A) = {u € Hl{Q.); cd^u £ H^{n)} . 

Let T > 0. We shall use the following notation i}' = U Q2, Q = 
(0,r) xn, Q' = (0,r) x n', Qi = (0,T) xQi, i = 1,2, r = {0,L}, and 
S = (0,T) X r. We also set S = {a'}. We consider the following parabolic 
problem 

idty + Ay = fm Q', 
\y{0,x) =yo{x) inn . 

(real valued coefficient and solution), for yo G L^(f]) and / € L?'{Q), with 
the following transmission conditions at a' 



(TC) 



(y{a'-)=y{a'+), 
\c(a'-)a,2/(a'-) = c(a'+)a,y(a'+). 



Now, we introduce finite-difference approximations of the operator A. 
Let = x'q < x'l < . . . < x'j^j^i = a' < . . . < x'^_^_.^j^i < x^_,_,^_,_2 = L. 
We refer to this discretization as to the primal mesh Tl := (x9i<i<n+m+i. 
We set := n + m + 1. We set h'._^i = x[_^_^ — x\ and x\_^y = {x'^j^^ + 

x9/2, i = 0, . . . , n + m + 1, and h' = maxo<j<„4-m+i ^' , i ■ We call 9Jt := 

*+ 2 

(^i+i)o<j<n+m+i the dual mesh and set h[ = x'^^i—x'.^^ = {h'.^i+h'.^i) /2, 
0, . . . , n + m + 1. 



In this paper, we shall address to some families of non uniform meshes, 
that will be precisely defined in Section 11.21 
We introduce the following notation 

[piAa = Pi(a+) - pi(a"), (1.3) 
Ma = P2(n+^)-p2(n + ^), (1.4) 

3 1 

[Pi-^p2]a = Pi(a+)p2(n + 2) ~ Pi(a~)P2(ra + 2)- (l-^) 

We follow some notation of [BHLlOa] for discrete functions in the one- 
dimensional case. We denote by and the sets of discrete functions 
defined on 9Jl and 9Jt respectively. If u € (resp. C^), we denote by Ui 
(resp. Uj_,_i) its value corresponding to x'^ (resp. 2;.^i). For u G we 

define 

n+m+l 

M _ 



u 



And for u G we define JqU ■= Jq u^{x)dx = Yl'i=i^^ h'-Ui. 
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For ueC^ we define 

n+m+l 

i=0 

As above, for u G , we define u := vF{x)dx = Er=o"^^ K+^'^i+h ' 
In particular we define the following L inner product on (resp. C^) 

{u,v)l2 = / u^{x)v^{x)dx, resp. {u,v)l2 = / vP^{x)v^{x)dx. 
Jn Jn 

For some u G C^, we shall need to associate boundary conditions = 
{uo,Un+m+2}- The set of such extended discrete functions is denoted by 
(^ajlua9Jt_ Homogeneous Dirichlet boundary conditions then consist in the 
choice uq = Un+m+2 = 0, in short = 0. We can define translation 

operators difference operator D and an averaging operator as the map 
^muam ^ ^-^^^ ^y 

{T~^u)i+i ■= Ui+i, {t~u)-_^_i := Ui, i = 0,.. .n + m + l, 

We also define, on the dual mesh, translation operators r^, a difference 
operator D and an averaging operator as the map — >■ given by 

{T'^u)i := u-_^_i, {T~u)i := u-_i, i = 1, . . . n + m + 1, 

{Du)i := ^(t^^ - T'u)i, u := ^(t+ + t~)u. 

1.2 Families of non-uniform meshes 

In this paper, we address non-uniform meshes that are obtained as the 
smooth image of an uniform grid. 

More precisely, let JIq =]0) 1[ and let i? : R — >■ M be an increasing map 
such that 

^{no) = n, i?gC°°, infi?'>Oand ^a) = a' (1.6) 

with a to be kept fixed in what follows and chosen such that a G (0, 1) fl Q, 
i.e a = I with p,q eW. Clearly, we have q > p. We impose the function -d 
to be affine on [a — S,a + 5] '&\[a-s,a+S] (for some 5 > 0). 

Given r G N* and set m = (q — p)r and n = pr. The parameter r is 
used to refine the mesh when increased. Set a = Xn+i = -Xpr+i- The interval 
^01 = [0, a] is then discretized with n = pr interior grid points (excluding 
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and a). The interval = [o-, 1] is discretized with m = [q — p)r exterior 
grid points (excluding a and 1). Let TIq = {ih)i<i<n+m+i with h = 
be uniform mesh of f^o and OKq be the associated dual mesh. We define a 
non- uniform mesh 9K of O as image of SUIq by the map i?, settings 



x'i = 'd{ih), Mi G {0, n} U {n + 2, n + m + 2} 
x'^+,:=a' = ^{a). (1.7) 

The dual mesh 9Jl and the general notation are then those of the previous 
section. 



1.3 Main results 

With the notation we have introduced, a consistent finite-difference ap- 
proximation of Au with homogeneous boundary condition is 



for u € C™^^^ satisfying 



A'^^u = -D{cdDu) 
u\dQ, = u^^^ = 0. We have 



{A'^u 



i+4 



1 Uj + l-Uj 

»+7 



, Ui—Ui-l 



hi 



1, .., n + m + 1. 



For a suitable weight function ip (to be defined below), the announced 
semi-discrete Carleman estimate for the operator = —dt + with a 
discontinuous diffusion coefficient c, for the non-uniform meshes we consider, 
is of the form 

2 
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L2{n) 



h 



-2 



2 

L2(f7) 



for properly chosen functions 6 = 9{t) and Lp = (p{x), for all r > To(T + r'^), 
< /i < /lo and r/i(ar)-i < eo, < a < T and for ah u G C°°(0,r;C^^) 
satisfying the discrete transmission conditions, where To,/io,eo only depend 
on the data. We refer to Theorem ST] (uniform mesh) and Theorem l5.6l (non 
uniform mesh) below for a precise result. The proof of this estimate will be 
first carried out for piecewise uniform meshes, and then adapted to the case 
of the non-uniform meshes we introduced in Section II. 2[ 
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Prom the semi-discrete Carleman estimate we obtain allows we deduce 
following weak observability estimate 

\q{0)\L^n) < Cobs lklli2((o,T)x^) \QiT)\l2(^n) > 

for any q solution to the adjoint system 

dtq + A'^q + aq = 0, q\on = 0. 

A precise statement is given in Section [6l 

Moreover, from the weak observability estimate given above we obtain a 
controllability result for the linear operator P'^. This result can be extended 
to classes of semi-linear equations 

{dt + A'^)y + Giy) = l^v, y £ {0,T) y\dn = 0, y{0) = yo, 

with G{x) = xg{x), where g € L°°(]R) and 

3 

\g{x)\ < KW'{e + \x\), x€R, with < r < -. 
We shall state controllability results with a control that satisfies 

\\v\\l2(q) <C\yo\ . 

Thanks to one space dimension the size of the control function is uniform 
with respect to the discretization parameter h. 

1.4 Sketch of proof of the Carleman estimate 

The main idea of the proof lays in the combination of the derivation of 
a discrete Carleman estimate as in [BHLlOal IBL12] and techniques used in 
[BDLOTj to achieve such estimates for operators with discontinuous coeffi- 
cients in the one-dimensional case. 

We set V = e-'^u yielding e'^Pe-'^v = e"^fi in Q'q if Pu = fi 

We obtain g = Av + Bv in Q'q, with A and iB 'essentially' selfadjoint. 

We write |j5||^2 = ||^v||^2 + |[i?f|[^2 + 2{Av, Bv)i2 and the main part 
of the proof is dedicated to computing the inner product {Av, Bv)i2(^q'^-^, 
involving (discrete) integration by parts. 

We proceed with these computations separately in each domain Qqi, 
fio2- As in [BL12j we obtain terms involving boundary points x = and 
x = 1 such as v{0),v{l),dtv{0),dtv{l), , (Dv)„_,.„+3 . In our case 

we obtain additional terms involving the jump point o such as v{a), dtv{a), 
Vn+llV^_^_3, {Dv)j^_^_i, {Dv)^_^_3. Main difficulties of our work come from 
dealing with these new terms. To reduce the number of terms to control, we 
find relations among connecting these various values at jump point allowing 
to focus our computations on terms only involving v{a),dtv{a) and {Dv)^,i. 



7 



Those relations are stated in Lemma 13.171 In the hmit h they give 
back the transmission conditions for the function v = e~^'^u used crucial 
way in |BDL07| . The idea of this technique comes from a similar technique 
shown in continuos case by |BDL07j . 

The discrete setting could allow computation on the whole il. Yet such 
computation would yield constant that would depend on discrete derivatives 
of the diffusions coefficient, yielding non-uniformity with respect to the dis- 
cretization parameter h. This explains why we resort to working on both 

and fli separately and deal with the interface terms that appear. As in 
[BDLOT] the weight function is chosen to obtain positive contributions for 
these terms. 

Sketch of proof Theorem 

1. We compute the inner product (Av, Bv) in a series of terms and collect 
them together in an estimate (see Lemma l4.4I ^Lemma I4.12p . In 
that estimate, we need to tackle two parts: volume integrals, integrals 
involving boundary points and the jump point. Volume integrals and 
boundary terms are dealt with similar to |BL12j . Terms at the jump 
point require special case. 

2. Treatment of terms the jump point 

• Terms at jump point involving dtv : when treating the term Y13 
we obtain a positive integral of {dtv{a))'^ in the LHS of the esti- 
mate as shown in Lemma 14.151 We keep this term in the LHS 
of the estimate. 

• Other terms: We collect together the terms at the jump point 
that already exist in the continuous case. As in [BDLOT] we ob- 
tain a quadratic form because of the choice of the weight function 
(jump of its slope). This allows us to obtain positive two inte- 
grals involving f^(a), {Dv)'^ i in the LHS of our estimate (see 

Lemma I4.14p . 

• The remaining terms at the jump point are placed in the RHS 
of estimate. After that, we apply Young's inequality to them (as 
shown in Lemma I4.16P and they then can be absorded by the 
positive integrals involving f^(a), {Dv)'^ ^ , {dtv{a))'^ in the LHS 

of estimate as described above. 
1.5 Outline 

In section [21 we construct the weight functions to be used in the Carle- 
man estimate. In section [3] we have gathered some preliminary discrete cal- 
culus results and we present how transmission conditions can be expressed 
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in the discretization scheme. Section [5] is devoted to the proof the semi- 
discrete parabohc Carleman estimate in the case of a discontinuous diffu- 
sion cofficient for piecewise uniform meshes in the one-dimensional case. To 
ease the reading, a large number of proofs of intermediate estimates have 
been provided in Appendix A. This result is then extended to non-uniform 
meshes in Section [5j Finally, in Section [Gj as consequences of the Carleman 
estimate, we present the weak observability estimate and associated some 
controllability results. 



2 Weight functions 

We shall first introduce a particular type of weight functions, which are 
constructed through the following lemma. 

We enlarge the open intervals 1^2 to large open sets Oi, ^2. 

Lemma 2.1. Let 0,2 be a smooth open and connected neighborhoods of 
intervals Qi, ^^2 ofM and let uj <zVt2 he a non-empty open set. Then, there 
exists a function tp G C{Cl) such that 



ipi in Qi, 
i>2 in VL2, 



with ipi e C°°{Qi),i = 1,2, V > m 0, -0 = on T, ■02 / in ll2 \ w, Tp[ / 
in and the function if) satisfies the following trace properties, for some 
ao > 0, 

{Au,u) > ao u € M^, 
with the matrix A defined by 

A = 

with 



an ai2 
0-21 022 



On = [V''*]a', 

022 = [ci^'Ml'{i^'){0^) + [c\4^'fMa', 
012 = 021 = [cip'-k]a'{lp'){o'^), 

(see the notation ()1.3p - (jl.Sp introduced in Section \l.l\) . 

Remark 2.2. Here we choose a weight function that yields an observation 
in the region lo <zVt2 in the Carleman estimate of Section^ This choice is 
of course arbitrary. 

Proof. We refer to Lemma 1.1 in [BDL07j for a similar proof. □ 
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Choosing a function tp, as in the previous lemma, for A > and K > 
1 1 I loo' define the following weight functions 

(^(x) = e^^(^) - e^^ < 0, ^{x) = e^^^''\ (2.1) 
r{t,x) = e"W^(^), pit,x) = {r{t,x))-\ 

with 

s{t)=Te{t), r>0, e{t) = {{t + a){T + a-t))-\ 

for < a < T. 
We have 

max6l = 61(0) = 9(T) = a'^iT + a)-\ (2.2) 

[0,T] 

and min^ > T^^. We note that 

[0,T] 

dtO = {2t - T)e^ . (2.3) 

For the Carleman estimate and the observation /control results we choose 
here to treat the case of an distributed-observation in w C 0,. The weight 
function is of the form r = e*"^ with cp = e'^'^ , with tp fulfilling the following 
assumption. Construction of such a weight function is classical (see e.g 
[Fl96]). 

Assumption 2.3. Let lo C Q be an open set. Let 0, be a smooth open and 
connected neighborhood of Cl in R. The function ip = 'ijj{x) is in (y(Q,R^, 
p sufficiently large, and satisfies, for some c> 0, 

ip > in 0, |V'i/'| > c in ^\uJo, 

dn^{x) < -c < 0, dltpix) < in Van. 

where Vqq, is a sufficiently small neighborhood of d^l in 0,, in which the 
outward unit normal n to Q is extended from dO,. 

3 Some preliminary discrete calculus results for 
uniform meshes 

Here, to prepare for Section [H we only consider constant-step discretiza- 
tions, i.e., h-,1 = h, i = 0, . . . ,n + m + 1. 

l-f- 2 

We use here the following notation: JIq = (0)1); ^oi = (0, a), ^02 = 
(a, !),% = Ooi U Qo = (0, T) X 1^0, Q'o = (0, T) x %, Qoi = (0, T) x J^oi 
with i = 1, 2 and 80.0 = {0, 1}. 

This section aims to provide calculus rules for discrete operators such 
as Di , Di and also to provide estimates for the successive applications of 
such operators on the weight functions. To avoid cumbersome notation we 
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introduce the following continuous difference and averaging operators on 
continuous functions. For a function / defined on we set 

r+/(x) := fix + h/2), t-/(x) := f{x - h/2), 
Bf{x) := (t+ - T-)fix)/h, fix) = (t+ + T-)/(x)/2. 

Remark 3.1. To iterate averaging symbols we shall sometimes write Af = 
f , and thus A^f = f. 

3.1 Discrete calculus formulae 

We present calculus results for finite-difference operators that were de- 
fined in the introductory section. Proofs can be found in Appendix of 
[BHLlOa] in the one-dimension case. 

Lemma 3.2. Let the functions fi and f2 be continuously defined in a neigh- 
borhood ofCl. We have: 

D(/l/2)=D(/i)/2 + /lD(/2). 

Note that the immediate translation of the proposition to discrete func- 
tions /i , /2 € C^^ and 91,92 e C®^ IS 

Dififi) = Difi)f2 + fiDif2), Digig^) = ^(51)52 + 51^(52). 

Lemma 3.3. Let the functions fi and f2 be continuously defined in a neigh- 
borhood ofQ. We have: 

/D'2 = /l/2 + ^D(/i)D(/2). 

Note that the immediate translation of the proposition to discrete func- 
tions /i , /2 € and IS 

112= fih + ^Dif,)Dif2), gm = 9-19-2 + ^^(51)^(52). 

Some of the following properties can be extended in such a manner to 
discrete functions. We shall not always write it explicitly. 
Averaging a function twice gives the following formula. 

Lemma 3.4. Let the function f be continuously defined over M. We then 
have 

A^f :=/=/+^DD/. 

The following proposition covers discrete integrations by parts and re- 
lated formula. 



11 



Proposition 3.5. Let f € i^p^^dWl g ^ (j-<OT_ have the following 

formulae: 

/ f{Dg) = - {Df)g + fn+m+2gn+m+l " foQ^, 
/ fg = I fg - -^fn+m+29n+m+^ ~ nfo9^- 

Lemma 3.6. Let f he a smooth function defined in a neighborhood of VL. 
We have 

T^f = f±^ r d,f{.±ah/2)da, 
^ Jo 

A^f = f + Cjh^ I ^{l- W\)dlfi. + ljah)da, 

D^f = dif + C'jh^ (1 - \a\y+^di:^^f{. + l,ah)da, j = 1, 2, h = \^h = 1- 

3.2 Calculus results related to the weight functions 

We now present some technical lemmata related to discrete operators 
performed on the Carleman weight functions that is of the form e^'^ , = 
gAi/j _ g-^^^ where ^ satisfies the properties listed in Section [2] in the domain 
Oq. For concision, we set r{t,x) = e^^^^'^^^^ and p = r^^, with s{t) = T6(t). 
Prom Section [21 we have fp\^^^ = ^iln^^' ^boi ~ ^^In^^ where ipi € C'^{(loi). 
Then p = e~^^ can be replaced by 

= e-'^f^ with If I = e^'^i - e^^ in domain Qqi 
p^ = e-"^^ with (^2 = e^^^^ _ ^XK (joj^ain 1^02 

And r = p'^^ is also replaced by 

ri = p^^ in domain JIqi 
^2 = i'^ domain Qq2 

The positive parameters r and h will be large and small respectively and we 
are particularly interested in the dependence on r, h and A in the following 
basic estimates in each domain i^oi, r2o2- 
We assume r > 1 and A > 1. 

Lemma 3.7. Let a, /3 G N, i=l,2. We have 

+ a/3(s(/.i)"A"+''-iOA(l) + s"-ia(a - 1)0a(1) = Oa(s"). 
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Let a € [—1, 1], we have 

(r,(t, .)id^p^)it, . + ah)) = Oa(s"(1 + ishf))e'>^^'^\ 

Provided < Th{max[Q 'j^j9) < R we have dx{ri{t, .){d'^ pi){t, . + ah)) = 
CA,.ft(s'"')- The same expressions hold with r and p interchanged and with s 
changed into -s. 

A proof is given in [BHLlOal proof of Lemma 3.7] in the time indepen- 
dent case. Additionally, we provide a result below to the time-dependent 
case whose proof is refered to |BL12l proof of Lemma 2.8]. Note that the 
condition < Th{max^Qrp^9) < ^ implies that s{t)h < ^ for all t € [0, T]. 

Lemma 3.8. Let a € N, i=l,2. We have 

dtirid^Pi) = s'^reoxii). 

With Leibniz formula we have the following estimates 
Corollary 3.9. Let a,/3,5 e N, i=l,2. We have 

diirfid^p,)dlpi) = (a + /3)^(-.^,)"+/^A"+^+^(vV'.)"+''+' 
+ 6{a + /3)(s,/.i)°+^A"+'^+'^-iO(l) 
+ s"+^-\a{a - 1) + /3(/3 - 1))Oa(1) = Oa(s"+^). 

The proofs of the following properties can be found in Appendix A of 
IBHLlOaj . 

Proposition 3.10. Let a € N, i=l,2. Provided < r/i(ma3;[o^T]^) ^ ^! we 
have 

r.T^d^Pi = nd^pi + s'^OxA^h) = s"Oa,j?(1), 
r.A^d^p, = rid^pi + s'^OxAsh)^ = s'^OxA'^), k = 0,l, 2, 
nA'^Dp, = r.d^pi + sOxA^h)^ = sOxA^)^ k = 0,l, 
nD^p, = r.d^^p, + s^OxA^hf = s^OxA^)- 

The same estimates hold with pi and rj interchanged. 

Lemma 3.11. Let a,/3 G N and k = 1,2; j = l,2;i = 1,2. Provided 
< T/i(max[o,T] ^) ^ ^> w^ have 

D\dP,{r,d^Pi)) = d'^+^ind^p,) + h^OxAn, 
A^d^ind'^Pi) = ^^{ri^'^p^) + h^OxAn- 

Let a € [-1, 1], we have D^d^{ri{t, .)d''pi(t, . + ah)) = OxA^^"^)- The 
same estimates hold with ri and pi interchanged. 
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Lemma 3.12. Let a, f3,5 e N and k = l,2;j = l,2;i = 1,2. Provided 
< Th{max[o^T]G) < A, we have 

D'^di{rf{d^Pi)d^,p,) = d',+\rHd^+^Pi) + h'0,,^{s^+^) = 0,,^{s^+^). 

Let a, a' G [—1,1]- We have 

AW'{r,{t, .fid'^piit, . + ah))d^p,iy, . + a'h)) = OxA'^""^^), 
D^d\n{t, .)\d^pi{t, . + ah))d^pi{t, . + a'h)) = OxM^"^^)- 

The same estimates hold with Vi and pi interchanged. 

Proposition 3.13. Let a G N and A; = 0, 1, 2; j = 0, 1, 2; i = 1, 2. Provided 
< sh < A, we have 

D^A^d'^inD^i) = d^+'^ind^pi) + sOx,A{shf = sOx,d^), 

D\nD''pi) = dlind-'pi) + s-'Ox,^{shf = s'OxM^), 

nA^Pi = 1 + Ox,^{sh)\ D^nA^pi) = Ox,^{sh)\ 
The same estimates hold with ri and pi interchanged. 

Proposition 3.14. Provided < r/i(maa;[o^T]^) ^ CLf^d a is hounded, we 
have 

dt{r,{.,x){dyi)[.,x + ah))=Ts'^e{t)Ox,^{l), 
dt{r,A^Pi) = T{sh)^e{t)OxAl), 
dtinD^Pi) = Ts^9{t)Ox,sia)- 

The same estimates hold with ri and pi interchanged. 

Proposition 3.15. Let a, f] e N and k = 0, l,2;j = 0, l,2;z = 1,2, 

provided < sh < ^, we have 

AW''df'irUd'')D^d = a^^(rf(5V)5ft) + s"+'OA,^(s/i)' = s"+iOA,^(l), 
AWdf'{rUd'')A''pi) = d',+P{n{d'^pi)) + s'^OxMi^h)^) = s'^OxA^), 

AW^d^irUd^D'^Pi) = 9^^(rf + s'^+^OxAsh)' = s"+'Oa,j^(1), 

and we have 

AW'^d^irfDpiD-'poi) = 5^"(r2(5ft)9V0 + s'OxA^h)'' = s^Ox^^), 
AW^cPirfDpiA^Pi) = d^.+'^indpi) + sOxM^h)-" = sOx,d^). 



14 



3.3 Transmission conditions 



We consider here discrete version of the transmission conditions (TC) at 
the point a. For u G we set / := D{cdDu) we then have 

u{a~) = u{a+) = Un+i, 

Remark 3.16. These transmission conditions provide the continuity for u 
and the discrete flux at the singular point of coefficient up to a consistent 
factor. 

Prom these conditions, we obtain the fohowing lemma whose proof is 
given in Appendix A 

Lemma 3.17. For the parameter A chosen sufficiently large and sh suffi- 
ciently small and with u = pv we have 

[i<CdDv\a = (cd£'w)„+3 - {cdDv)^^)_ = JlVn+1 + J2icdDv)„^l + Js.h{rf)n+i (3.1) 

where 

Ji = (1 + OxAsh))H*cH']a + sOx.siish), 

J2 = Ox^^ish), J3 = (1 + OxAsh)). 

Furthermore, we have 

dtJi = sTe{t)OxAsh). 

dtJ2 = Te{t)OxAsh), dtJs = Te{t)OxAsh). 

For simplicity, (jS.ip can he written in form 

[kCdDv\a = \s[^ccj)ii']aVn+i + ro, (3.2) 
where tq = XsOx,fi{sh)vn+i + Ox,ii{sh){cdDv),^^i + hil + Ox,<;i{sh)]{rf)n+i. 



4 Carleman estimate for uniform meshes 

In this section, we prove a Carleman estimate in case of picewise uniform 
meshes, i.e, constant-step discretizations in each subinterval (0, a) and (a, 1). 
The case of non-uniform meshes is treated in Section [5l 

We let ujQ C ilo2 be a nonempty open subset. We set the operator to 
be V'^ = -dt + A'^^ = -dt - D{cdD), continuous in the variable t £ (0,r) 
with r > 0, and discrete in the variable x (z Qq. 

The Carleman weight function is of the form r = e^^ with ip = e^^ — e^^ 
where ij: satisfies the properties listed in Section [5] in the domain VIq. Here, 
to treat the semi-discrete case, we use the enlarged neighborhoods fioii ^02 
of J7oi) ^^02 as introduced in Lemma [2Tl This allows one to apply multiple 
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discrete operators such as D and A on the weight functions. In particular, 
we take ip such that dxip > in Vq and dx'il^ < in Vi where Vq and Vi are 
neighborhoods of and 1 respectively. This then yields on OVIq 



{rDp)o < 0, irDp)n+m+2 > 



(4.1) 



Theorem 4.1. Let ivq C r2o2 be a non-empty open set and we set f := 
D{C(iDu). For the parameter A > 1 sufficiently large, there exists C, tq>\, 
ho > 0, eo > 0, depending on ojq so that the following estimate holds 

2 







2 

+ r 










LHQo) 









+h- 



i2e'"^u 



L^{{0,T)xuo) 



U\t=0 



+ h~ 



U\t=T 



2 



(4.2) 



for all T > tq{T + T^), < h < ho and Th{aT) ^ < eo and for all u G 
C7°°(0,r;C^) satisfying u\Qa^ = 0. 

Remark 4.2. Observation was chosen in ilo2 here. This is an arbitrary 
choice (see Remark \2.S\) . 

Proof. We set /i := -P™ = dtu + D{cdDu) and / = D{cdDu). At first, 
we shall work with the function v = ru, i.e., u = pv, that satisfies 



We have 



r(dt{pv) + D{cdD{pv))j = rf, in Q'^. 
rdt{pv) = dtv + r{dtp)v = dtv — T{dt9)tpv. 



(4.3) 



We write: g = Av + Bv, 

where Av = Aiv + A2V + ^311, Bv = Biv + B2V + B^v with 

Aiv = rpD{cdDv), A2V = cr{DDp)v, A^v = —T{dt9)ipv, 

Biv = 2crDp Dv, B2V = —2sc(j)"v, B^v = dtv, 

h _ _ 

9 = rfi- -rDp{Dcd){T+Dv - Dv) - —{Dcd)r{DDp)Dv 

-hOiiylTp - {r{Dcd) ITp + hO{l)r{DDp))v - 2 

as derived in [BL12j . 

Equation ()4.3p now reads Av + Bv = g and we write 



sc 



\\M?mQ[;) + \\Bv\\12(^q,^^ + 2{Av,Bv)l2^q,^) = \\g\\l2^Q,^ 



(4.4) 



First we need an estimation of ||5||£^2(q/ ). The proof can be adapted from 
[BHLlOa] . 
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Lemma 4.3. For T/i(max[o,T] d) < ^ we have 

Ml^iQ'^) < CxMWrfiWhiQ',)) + P^lli^wj.) + \\sDv\\l,^Q,^^ . (4.5) 

Most of the remaining of the proof will be dedicated to computing the 
inner product {Av, Bv)i2(^Q'^^y Developing the inner-product {Av, Bv)]^2^qi^-^, 
we set lij = {AiV, Bjv)i2(^Q/^y The proofs of the following lemmata are 
provided in Appendix A. 

Lemma 4.4 (Estimate of In). For Th{max^Q^T]G) < R we have 

hi>- ! sX'{c''cl>{^'fUDvf - Xn + Yn, 
Jq' 

where Xn = Jq, i'ii{Dv)'^ with un of the form sX(f)0{l) + sO\^^{sh) and 
ill — -I- -111 ^ '^ii ' 

T 

2 

3 



yS^ = {l + OxAsh)){cca){l){rD-p){l){Dv)] 
(1 + Ox,^{sh)){ccd)mr'D-pm{Dv)\, 
y}i'^ = £ sXct>{a)-c,{a)[{c^'){a+)iDv)l^, - {c^'Ka-){Dv)l^,) , 
^11 ''^= f sOxAsh)\Dv)l^,- r sOxAsh?{Dv) 

Jo 2 Jq 



^2 

3 ■ 



Lemma 4.5 (Estimate of Ii2). For T/i(maa;[o,r]^) ^ -^^ ihe term I12 is of 
the following form 



112 

'Q'o 



2 [ sX\c^(t>iil^'fUDvf - X12 + Y12, 



with 



Y12 = [ sX^^{a)v{a)[cii;'f*CdDv]a 
Jo 

+ 5i2v{a){cDv)^^^+ 8i2v{a){cDv)^^i_, 
where (512,(512 are of the form s{X(f){a)0{l) + 0\^^{shY) and 
X12 = / vi2{Dvf + [ sOxA^)vDv, 

Jq'o Jq'o 

where 

z/12 = s\(i>0{l) + sOxAh + {shf). 
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Lemma 4.6 (Estimate of /13). There exists ei(A) > such that, for < 
r/i(max[o,T]^) ^ ^iW) the term I13 can be estimated from below in following 
way: 

hz>- I ex,^{l){Dv{T)f - Xis + Ha- 



with 



Xi3= I {s{sh) + T{shfe)Ox,^{l){Dvf+f s-^Ox,^{sh){dtv) 

jQ'n JQo 



"0 

r-T rT 



Yi3 = - rp{a'^)dtv{a)(cdDv)^j^^+l rp{a )dtv{a){cdDv)^j^i. 
Jo Jo 

Lemma 4.7 (Estimate of I21). For Th{max^Q t]6) < ^, the term I21 can be 
estimated as 



'Q'o 

with 



/2I > 3 / AS303c2(^')4y2 _ ^ Y21, 

JQU JO'r. 



where 

/X21 = {s\<t>fO{l) + s''Ox,^{l) + s^Ox,A{shf, U2i = sOx,A{sh)\ 

and 

121 — 121 + ^21 + -^21 + -^21 ' 

Jo 2 Jq 2 

Jo 

Jo 2 2 

Jo 

Lemma 4.8 (Estimate of 122)- For sh < A, we have 

I22 = -2 / c2s3aV^(V'')V - X22 + I22, 

JQ'r, 
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with 



^22 — ^22 + -'^22 ' 



and 



Jo 

X22 = / /^22'y^ + / J^22{Dvf, 
JO' JO.' 



where 

/X22 = {sX^fOil) + S-'O^^^il) + 530a,^(s/i)2, i.22 = sO^,^{shf. 

Lemma 4.9 (Estimate of J23). For r/t(max[o_r]^) ^ the term I23 can be 
estimated from below in the following way 

I23 > I ^ s\Ox,^{l)vl^^ + Ox,a{1)vI^^) - X23 + Y23, 

*f fin 



with 



X23 = [ Ts^eOxA^y + [ s-^Ox,^{shf{dtvf 
JQo JQo 

+ [ {shfsOxA'^){Dvf, 
JO'. 



and 



^23 - J:^23 + ^23 + ^23 ' 



Jo 

Yi^ = 0,,^{sh)'v\a)\lzl 
Lemma 4.10 (Estimate of /31). For Th{max^Q^x]G) ^ ^> we have 

I31 = —X31 + Y31 , 



2 



with 

X31 = [ Tes^Ox,^{iy + / TeoxA^hf{Dvf, 

JQ'o JQ'o 
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and 



^31 - + ^31 ' 



Yg^= r Tes^Ox,^{l)v\a). 
Jo 

Lemma 4.11 (Estimate of /32)- IB LI 2^ For Th{max[QT]0) < ^, the term 
I32 can be estimated from below in the following way 



'32 



Lemma 4.12 (Estimate of 133). WLISI proof of Lemma 3. 9j For Th{max[QX]G) ^ 
^, the term L33 can be estimated from below in the following way 



^ JQ' 



'Q'o 

Continuation of the proof of Theorem 4.1. Collecting the terms we have 
obtained in the previous lemmata, from ()4.4p and ()4.5p for < rh^max^o x]^) < 
ei(A) we find 

2/„2^^„/,/^l2^^n„,^2 I o / „2„3 \4^3/„;,/^4^2 



< 



\\MrmQ'„) + ||i?«|li.(Q.) + 2 / s\\c'<p{^/YUDvY + 2 / c^s^AV(^') 

''Q'o Q'o 
+2(nW + 4^'^)) + 2(yA'-'^ + 4'^''^) + 2^13 

C^A,^(ll'-/illi^(Q^)+ /,^'(Co+CJ+ j {Dv{T)f)+2X + 2Y 

J On ^ Or, 



with 



Y=- (y/i'^'' + Y12 + Y^l^^^^ + + Y22 + 1^23 + ^31 

X = Xii + Xi2 + Xi3 + X21 + ^22 + ^23 + ^31 + ^32 + ^33 

+CxA\\^v\\l.iQ'o)+h'\\sDv\\l.^Q,^). 

With the following lemma, we may in fact ignore the term Y^^ + Y^l''^^ in the 
previous inequality. 

Lemma 4.13. For all X there existsO < e2(A) < ei(A) such thatforO < Th{max\i^T](^) ^ 
£2 (A) we have ' + Y^l'^'^ > 0. 
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Recalling that 7-0 > C > in il\ujQ we may thus write 



II^Hli2(Q,) + WBvWi.^Q,^ + / s{Dvf + / 

•I Q'q Q'o 



'J UJQ '^0 J UJQ 



+ [ ''H-o+^U-t)+ f {Dv{T))A+2X + 2Y. (4.6) 
Jn' Jn' / 



Lemma 4.14. With the function ijj satisfing the properties of Lemma \2.1\ and for 
Th{maX[o^T](^) < ^, we have 



Yii^^ + '"'^ > Can / s\(t){a){cdDv)l + Ca^ / s"" f {a)vl+^ + Mi + M 



"+5 

'Jo 



T 



with ao as given in Lemma \2.1\ and where 

Mr = / sOx[l)rl+ j s^Ox{l)roVn+i+ [ sOx{l)ro{cdDv)^+i 
Jo Jo Jo ^ 

+ / s^CA,j?(s/i)w^+i + / sOA,^(s^)w„+i(c<ji:)u)„_^i + / sOx,si{sh)roVn+i, 
Jo Jo Jo 

with To as given in Lemma \3.17\ and 

Ml — Ml + Ml I 

w/iere 

Jo ^ Jo ^ 

Mp'' = / Ox^si{sh)(cdDv) ^j^3Vn+i + / s^Ox,si{sh){cdDv)^j^iVn+i- 
Jo Jo 

For a proof see Appendix A. 
Lemma 4.15. With < e3(A) < e2(A) sufficiently small we obtain 

Yi3> [ Cx,sih{dtv{a)f + [ {sTeOxMsh)+T^0^OxAsh))v^a) 
Jo Jo 

+ sOx.si{lV{o)\\Zo + I Oxs.{sh)dtv{a){cdDv)^^i + [ OxAmv{a)h{rfi)n+i. 

Jo ^ Jo 

where Cx,si is positive constant whose value depends on A and sh. 
For a proof see Appendix A. 

If we choose A2 > Ai sufficiently large, then for A = A2 (fixed for the rest of the 
proof) and < Tft,(maxjo,T] 6*) ^ £3, from (14.6^ and Lemma [4.141 and Lemma [4.151 
we can thus achieve the following inequality 
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\\M\l^Q'„) + \\Bv\\l.^^Q,^)+ f s\Dv\'dt+ f s\'dt 

J Q'o J Q'o 

+ Cao [ sicdDv)l^,+Cao [ sV(a) + / Cx,fMdtv{a)y 
Jo ^ Jo Jo 



S V 

^0 Jloq Jloq 



+ -'H-o {Dv{T)r + sv\a)\lzl 
J% Jn'„ ^ 

T pT 
2/)2/ 







+ f Ox,fi{l)dtv{a)h{r,h)„+i+2X + 2Y + 2Z , (4.7) 

Jo 

where Z_ ^ fir + fJ-i with fi^ and fii are given as in Lemma 14.141 and where 



X = I P-v^ + v{Dvf 



Q'o 







+ Xi2 + Xi3 + X23, + X31 + ^32 + -'^^33 1 

with /x = s'^0\^^{l) + s^O\,si{sh) and u of the form sO\^^{sh). 

By using the Young's inequahty, we estimate in turn all the terms of Y_, Z_ and 
the two terms at the RHS of (|4.7I) through the following Lemma whose proof can 
be found in Appendix A 

Lemma 4.16. For sh < ^, we have 



Jo Jo 

OxAsh)dtv{a){cdDv)^^+. <e [ OxA'^)HdMa))^+a I sOxAsh){cdDv)l^^. 







T T T T 







711 = sOxAshf Vii = OxAshf. 



\Yi2\< ai2<+i+/ Pi2h{dtv)i+^+ lAcdDv)f^^,+ rn2/i(r/i)^+i, 
Jo Jo Jo Jo 
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712 = OxA'^) V12 = OxAsh). 



(1,22) 



Jo Jo Jo Jo 

a2i = (s^OxA^) + sTH^OxAshf) P21 = OxA^hf, 
721 = sOxAshf V21 = OxAshf. 



-'22 



< 



^23 



< 



/ a22vl+i+ [ l322h{dtV)l^i+ [ j22{CdDv)l 1 + [ 7722/l(r/l)^+i, 

Jo Jo Jo ^ Jo 

a22 = [s^OxAshf + sT^e^OxA^hf) P22 - OxA^hf, 
722 = OxAshf V22 = OxAshf- 

I a23vl+,+ [ P2zh{dtv{a)f+ f ^2s{cdDv)l ,+ f n22h{rh)l+^, 
Jo Jo Jo Jo 

a23 = (s'OxA^) + sT^Q-'OxAshf) P23 = OxA^h), 
723 = sOxA^hf r?22 = OxA^hf. 



< a3ivl+i+ I323h{dtv)l+i + l2z{cdDvf i+ I 7722/i(r/i)^+i, 
Jo Jo Jo ^ Jo 



T 

(1) 



asi = (s'^reOxAsh) + sT'^e'^OxAsh)) /323 = OxAsh), 
723 = sOxA^h) V22 = OxA^f^)- 

Hi< ( aiu^+1 + / ^rh{dtv)lj^^ + / 7i(cd-Dw)^ , 1 + / mK^h) 
Jo Jo Jo ^ Jo 

ai = (s'OxAsh) + sT^'^OxAsh?) A - OxA^hf, 
71 = sOxAsh) Vi = OxAshf. 

Hr< ( arvl+i + / 0rh{dtv)l^-^ + [ jricdDvf i + [ r]rh{rfi) 
Jo Jo Jo ^ Jo 



2 

n+l! 



2 

n+1' 
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-ir = (sOxAsh) + esC'A,j^(l)) Vr = 0,.x,si{sh). 
Futhermore, we can estimate the term in X12 as follows 

/ sOxA^)vDv < f sOxA'^)ii?+ I sOx,simDvf 

Jq'o Jq'o •'Q'o 

sOxA^)V\ + I sOxAWvf 



Q'o 



Q'o 



Qo 



sOxA^y + / sOxA^)iDv)^. 



Q'o 
2 



2/i3 



by Lemma O and as J^,^OxA^) \^\ = Jug'^xA^)^ 
Observe that 

1 < T^e and 19(^611 < CT 

We can now choose £4 and ho sufficiently small, with < €4 < e3(A2), < ho < 
hi{X2), and T2 > 1 sufficiently large, such that for r > r2(T + T^), < h < ho, and 
Th{maX[o^T]S) < £4, from (j4.7l) and Lemma 14.161 we get 



s\Dv\^ 



+ Cao I s{cdDv)^^i + Cao 



„3„,2 



1^,1+1 + Cx^si / hidMa)) 



< CxM[\\rh\\l^Q,)+ f I s{Dvf+[ f 

y J J ojQ Jo J ujQ 

+h-'{ f Co+ / CT)+*^'(°)lt=n + f OxAmrfi)l+i 



sOxA^y+ s-^OxAsh){dtvf + s'TOOxA^)^ 



Qo 



(4.8) 



where we used that [Dv)^ < Cft-^^((T+?;)^ + (r^w)^) and the last three terms whose 
integral taken on domain Qo come from the term in X12, X13 and X23 respectively. 
As T > T2(T + r^) then s > T2 > and furthermore we observe that 



s ^dtv 



LHQ'o) 



< C 



1 


2 


1 


2 


1 




s 2 Bv 


+ 

LHQ'o) 




+ 

LHQ'o) 


Dv 


LHQ'o) ^ 



< Cx^r.A \\BV 



3 


2 








+ 


s^Dv 


' ) 




L^Q'o) 




LHQ'o)'' 



We then add the following terms hs^v'^^i and Jo hs ^{dtv{a))'^ on both 



the right hand side and the left hand side of (|4.8p . This allows us to change the 
domain of integration from Qq to Qo for the discrete integrals on the primal mesh. 
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No additional term is required for discrete integrals on the dual mesh. For sh 
sufficiently small and s > 1 sufficiently large, these terms at the right hand side are 
then absorbed by the terms at the left hand side. More precisely, with < eo < 64 
sufficiently small and for t > T2(T + T^), < h < ho, and < Th{max\fi T](^) < eq 
we thus obtain 



s ^dtv 



2 



+ 



s{Dv) 



Qo 

2 



Qo 



< C^A,.U ||r/i||^.(Q^)+ / / s{Dv) 

\ Jo Juja 



Ju:o 



-h-'if <J+.0.,^(iy(a)|*=j). (4.9) 

J no J no / 



Now we shall estimate the term sOx^ix{l)v'^(a)\^^^. We have 

n+m+l 



\t=T \\L^(n, 



It follows that, as sh is bounded 



, = > hv'^i > h IIvl 

q) j\t = T — II lt = 



T|lL~(no) ■ 



T\\L^no) 



Similarly, we treat the term sO\,ii{l)v'^ {a)^^^^ as 



\sOx,si{^)v'^ia)u=o \ < C'a.j? h ^ \\v\,^o\\mno) 
Therefore, (|4.9p can be written as 

2 



1 




2 


3 


+ 




+ 




i^(Qo) 




i^(Qo) 





s^Dv 



L2((o,T)XLJo) 



2 



L2((o,T)xwo 
,,2 



|t = / |t = T; 

<j no 



We next remove the volume norm s 2 Dv 
proof of Theorem 4.1 in |BHL10a] we thus write 



by proceeding as in the 







2 

+ T 




2 






LHQo) 









2 

i=(Qo) 



L2((o,T)xt^o) 
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As we have max^ < t^-, we see that a sufficient condition for Th(nia.x9) < en 

[0,T] - ^" ^[0,T] ' - 

then becomes Th(Ta)~^ < eq- To finish the proof, we need to express all the terms 
in the estimate above in terms of the original function u. We can proceed exactly 
as in the end of proof of Theorem 4.1 in |BHL10aj . 



5 Carleman estimates for regular non uniform meshes 

In this section we focus on extending the above result to the class of non piece- 
wise uniform meshes introduced in Section 11.21 We choose a function satisfy- 
ing p.6p and further ^\[a-s,a+s] is chosen affine (for some S > to remain fixed in 
the sequel). The way we proceed here is similar to what is done in |,BHL10a| . In 
this framework, we shall prove a non-uniform Carleman estimate for the parabolic 
operator = —dt + on the mesh 9Jl by using the result on uniform meshes 
of Section m 

By using first-order Taylor formulae we obtain the following result. 
Lemma 5.1. Let us define C, G and C, G as follows 



h 



C+i = -^^i e {0,...,n + m + l} G = e {!,..., n + r7i+ 1} 



These two discrete functions are connected to the geometry of the primal and dual 
meshes DJl and 971 and we have 

< inf I?' < C+i < supi?', Vi e 0, . . . ,n + m + 1 



no 



< inf t?' < Ci < sup Vi e 1, . . . , n + TO + 1 



We introduce some notation. To any u £ (j^OTuaon^ associate the discrete 
function denoted by Q^°u G c^i^ouaOTo (^ggj^gfi on the uniform mesh SJIq which 
takes the same values as u at the corresponding nodes. More precisely, if m = 
i=i i[x' ^,x' jw,, we let 

n-t-m-t-l 
i=l 

and {Q^l'u)o = uq, {Q^^°u)n+m+2 = u„+„i+2- Similarly, for u e C^, u = 
1=1 Mx'^.x'^^,]u,+ i,weset 



n+m+1 

i=0 



The operators Q^j" and Q^" are invertible and we denote by Q^o ^^"^ their 
respective inverses. We give commutation properties between these operators and 
discrete-difference operators through the following Lemmata whose proofs can be 
found in ,BHL10a) . 
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Lemma 5.2. WHLlOal see the proof of Lemma 5.2] 

1. For any u £ C^'^^^^ and any v G , we have 

D{QZ"u) = Q^iCDul DQ^v = Q^«(C Dv) 

2. For any u £ fmyg 



Lemma 5.3. \BHL10a[ see proof of Lemma 5.3] 
For any u £ and any v £ CP'^ , we have 



szo 

Q'\-l L.|2 



On 



2 



Futhermore, the same inequalities hold by replacing by uj and flo by ujq, respec- 
tively. 

For any continuous function / defined on (resp. on ^o) we denote by 
Hot/ = if{xr))o<^<n+m+2 € C^^OOT ^he sampling of / on (resp. U^J = 
ifiih))o<^<n+m+2 € C^oUOOTo the Sampling of / on DJlo). 

Lemma 5.4. ]BHL10a\ see the proof of Lemma 5.4] 
Let f be a continuous function defined on Q 

Q^;°noji/ = na,i„(/o^). 

In particular, for u £ C'-*^^^^ we have 

QZ'({T\^if)u) ^n^,{f o^){Q^>u). 

Moreover, by making use of Taylor formulae we get the following result 
Lemma 5.5. With Q defined as in Lemma \5.1\ we have 

\\DDi^\\ < oo, \\Di>\\ < oo, < llz^n , ||i>|l < oo 

II lloo ' II lloO ' IlllOO'llllCXj 

where v := -w^- 

Proof. From the definition of C, 2^° ^^'^ ^ acting on C'"''", D acting on C™" we 
have 



(5.1) 
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We find 



X, 1 = 



i~l 



K+1 = 



K-1 



^{{i + 2)h)-d{lh) ^ ~ ^^ 

2 2 
d{ih) - - 2)h) _ - l?,_2 



2 2 
By using Taylor formulae we write 



"9i+2 = 


d, + (2/i)i9'^ + 


{2hf {2hf 


(2M\.(4) 


2 ' 6 ' 


24 * 




- {2K)i)'^ + 


{2hf {2hf 


(2M%(4) 




2 6 ' 


24 ' 


= 






+ 0{h^), 






' 6 ' 24 ^ 


+ Oih^). 



0{h% 
0{h% 



Thus we have 
h' . 



2/i^: + ^C + 0(/^'), 



K+i = '^hd[ 

hLi = 2/i^: 
From (|5.ip we obtain 



(2^ 

2 ■ 6 " ' 24 
2 6 24 



- , 1 , N 



where 



and 



Thus, we have 



< (infi?')" < oo, 



which proves the first result. Next, we proceed with the second result in the same 
manner as above. We have 



2h 



1 

2h^h 



'i-l 
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By using the computations of h'^_i, above we find 



-{2h)H" + 0{h^) 



< 



\\^"\L 

(inf 



< oo, 



which yields the second result. 

Moreover, with the properties of C shown as in Lemma l5. II we can assert 



0< ll^lloo:ll'^lloo<0«■ 



□ 



From Lemmata 15.21 - 15.41 we thus obtain the following discrete Carleman esti- 
mate for the operator = —dt — D{cdD.) on the mesh 

Theorem 5.6. Let u <Z ^l2 be a non-empty open set and we set f :— D{cdDu). 
For the parameter A > 1 sufficiently large, there exists C, Tq > 1, ho > 0, eo > 0, 
depending on tu such that for any mesh 371 obtained from by (jl.6p - (jl.7p . we 
have 



< C 



LHQ) 

2 



2 



I 1-2 I TOip I r 

+h |e '^?^|t=o|i2(fj) 



L2((o,T)xw) 



2 



, (5.2) 



for allT > tq{T+T^), < h < ho andTh{aTy^ < co and for allue C°°(0, T; C^^^) 
satisfying — 0. 



Proof. We set w — Q^i°u defined on the uniform mesh OTq. By using Lemma [5.21 
we have 



Dw 



(5.3) 



We observe that the right-hand side of (|5.3I) is a semi-discrete parabolic operator 
of the form P^" = Ci-dt - jrDi^dD.)) applied to w, where 



We set := 77 = .iL _ and we find 



9 = iy + h^DDv + h^O{l), 



(5.4) 



by using Lemma 13.31 and Lemma 15.51 

Thus, the operator P^" can be written in form as 

P^^°w ^^'{-dtw- PDi^dDw) + h^O{l)D{£,dDw)). 

Moreover, using Lemma l3.2l we have 



vD{^dDw) = D{V^dDw) - D{v)^dDw. 
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We set P^°w := -dtw - D{D£^dDw) 
the properties of v and it follows that 

< hram <b < bmax and 



'dtW — D{bdDw) with — i>^d- From 



\D{bd)\ 



< oo. 



First, we shall obtain a Carleman estimate for Pq Then we shall deduce a 
Carleman estimate for the operator 



(5.5) 



Now, we consider the function i/i o : (t, x) M- ip(t,'d{x)Y By using the prop- 
erties listed in Lemma 12.11 and (|1.6p , we shall see that ip o i!) is sl suitable weight 
function associated to the control domain ujq — d~^{w) in J7oj i-e., that ^po■|) satisfies 
Lemma l2.1l for the domaims and uq- 

The important property to checking is the trace property. The remaining prop- 
erties are left to the reader. We set 



B 



hi bi2 

^21 &22 



with 



bn = [(^ot?)'*], 

622 = H-^ o o ^)'(a+) + [b^ib o d) 

&12 = &2i = [6(Voi?)'*]a(V°^9)'(a+) 



'3^ 



where b 



n^io c 



Morever, we have 'd'_^{a) = i9'_(a) and C„+i = K 



n+l 



(recall that 'd\[a-5,a+5] is an affine function). It follows that 



bii 

&22 



[^V'(^)^'*]^V''(^)(a+)^'(«+) + [{'-^ni^'i^WfMa 



&21 



[^V'(^?)^'*]aV''(^?(a+))^?'(a+) 



We can see that {Bw,w) ~ {Aw,w) > ao . This means that ipo-d satisfies 
the trace property. 

Through Theorem 14. 1[ we obtained a discrete uniform Carleman estimate for 



P^" and the Carleman weight function is of the form ro = 6^^^° , with ipo = ipo-ff 



e^^ where ipo = ip o -d on the uniform mesh DJIq . We can deduce the same 



result on the non-uniform mesh 971. Namely, through (14.21) we see that the following 
estimate holds 



< C 



2 


f r 


LHQo) 





i2e' 



'Dw 



2 

L2((o,T)xtuo) 



w 



2 



2 



(5.6) 
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and the constant C is uniform in h for r sufficiently large and with Th(aT)~^ < cq, 
for eo sufficiently small. Note that, setting cq — (inf^o i9')^0i we see that the 
condition rh' {aT)~^ < Iq on the size of the non- uniform mesh 3Jl implies the 
condition Th{aT)~^ < eg for the uniform mesh dJlo. 

From (|5.5p - (j5.6l) we deduce the following Carleman estimate for P^" 



0- 



2 

+ r 




2 


f r3 




2 


i^(Qo) 




i^(Qo) 









< (7(||gre^opOToy;| 



i"(Qo) 



L2((o,T)XLdo) 



2 



(5.7) 



Now, by using Lemma 15.51 we estimate ||e^^'^''£)(f>)^^£)y;||^^^^ ^ jj^ the RHS 



of (|5?71) as 



We see that 



Hence we find 



< C 

< C 



< C 

< C 



2 



< C 

< 

by using (|5.4p and Lemmata 15.21 - 15.4 



LHQo) 
2 



We treat \\e''fif+{cdDu)\\^,^Q^ 
similarly). We find 



(the term ||e*'^f_ (crfDu) 



can be treated 



\\e''^f+{cdDu)\\^,^Q^ - \\rf+{cdDu)\\^,^Q) < \\{T^r){cdDu)\\^,^Q^ < C \\{T^r)Du\\^,^Q) 

(5.8) 

We have T_r — r{pT^r) = r{l + 0\,^{sK)) (due to Proposition 13.101) . From 
that we can write 



which allows one to absorb by the term at the LHS of the Carleman estimate by 
choosing r sufficiently large. 
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Next, we estimate h'^ \\e'^'^'^o D{^dDw) 

|2 



in the RHS of (O) as 



LHQo) 



< Ch^ 

< Ch^ 

< Ch^ 

< Ch^ 

< Ch^ 

< Ch 



LHQo) 
2 

LHQo) 

2 

+ 

2 



LHQo] 
2 



2 



LHQo) 



by using (j5.4p and Lemmata 15.21 - lOl We proceed with an estimate as in (jS.Sp . 
We thus obtain 

which allows one to absorb by the term in the LHS of Carlenian estimate by choosing 
T sufficiently large. 

Futhermore, by using the previous Lemmata 15.11 - 15.41 and considering each 
term in (|5.7p separately, we see that we have the following estimates 



• For the first term in LHS of ((5?71) 





2 




LHQo) 



m 



> (sup^?') 



LHQo) 



LHQ) 



and a similar inequality holds for 



2 



• For the second term of LHS of (|5.7I) we use Lemma 15.21 and Lemma 
follows 



as 



LHQo) 



2 

LHQo) 



an, 



LHQo 



> 



• By using (|5.3p and Lemma 15.31 we have 



lL2(Qo) 



L^(Qo) 



LHQo) 
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• For the third term of RHS of ([577| 

|2 



e t=o 



L2(n) 



and a similar inequahty holds for e'^^'^°w|t=T 



• Finally, since i?(a;o) = uj we have 

2 

L2((0.T)xtJo) 



< 



The proof is complete. 



L2((o,T)xa;o) 



L2((o,T)xwo) 



L2((o,T)xw 



□ 



6 Controllability results 

The Carleman estimate proved in the previous Section allows to give observ- 
ability estimate that yields results of controllability to the trajectories for classes 
of semi-linear heat equations. 

6.1 The linear case 

We consider the following semi-discrete parabolic problem with potential 

dty + A^'^y + ay^l^v, t e {0,T) y\on = (6.1) 

The adjoint system associated with the controlled system with potential (|6.ip 
is given by 

-dtq + A^^'y + ay^O, t e {0,T) q\an = (6.2) 

We assume that a piecewise diffusion coefficient c satisfies (11.21) and 57 — (0, 1). 
From Carleman estimate (j4.2l) we obtain a following observability estimate. 

Proposition 6.1. There exists positive constants Co,Ci and C2 such that for all 
T > and all potential fucntion a, under the condition h < min(ft,o, ^1) with 



hi ^ Coil 



1 
T 



Mir 



any solution of (|6.2p satisfies 

l9(0)lL2(n) 



|9(T)I 



(6.3) 



with Cobs 



,C2(l+4-+T||a||^ + ||a|||, 



Remark 6.2. In comparision the observability inequality in continuous case which 
performed in \BDL07^ , we find that the observability inequality obtained here is weak 
since there is an additional term depending upon h at right-hand- side of inequal- 
ity (lOl) . 
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From the result of Proposition 6.1 we deduce the fohowing controUabihty result 
for system (|6.ip . 



Proposition 6.3. There exists positive constants Ci, C2, C3 and for T > a map 
LT,a ■■ R^' ^ L2(0,T;M®^) such that if h < min(/io,/i2) with 

hi^Co{l + ^+T\\a\\^ + \\a\\iy' 

for all initial data yo G M'^^ there exists a semi-discrete control function v given by 
V — La{yo) such that the solution to (|6.ip satisfies 

|2/(r)|^.(,,)<Coe-^^/''|2/oL2(^) 

and 

||w|Il2(q) < Co lyolL2(n) 

Note that the final state is of size e"'^^'' \yo\L^(^n)- The proof of these proposition 
are given in |BL12] . 

6.2 The semilinear case 

We consider the following semilinear semi-discrete control problem 

{dt+A''')y + Giy) = l^v, ye(0,r) y\an^O, y(0) = 2/0 (6.4) 

where w C il. The function G : M — > M is assumed of the form 

G{x)=xg{x), xeR, (6.5) 

with g Lipschitz continuous. Here, we consider the function g in two cases: g 6 
L°°(R) and the more general case as 

3 

|5(a;)| < inn''(e + |x|), x £ R, with 0<r<- (6.6) 

The results of semi-discrete parabolic with potential above allows one to obtain 
controllability results for parabolic equation with semi-linear terms whose proofs 
are given in |BL12| 

Theorem 6.4. We assume that g € L°°(M) and c satisfies (II. 2p . There exists 
positive constants Co,Ci such that for all T > and h chosen sufficiently small, 
for all initial data yo G MP^ , there exists a semi-discrete control function v with 

\Hl^Q) ^ C'lyolL2(o) 

such that the solution to the semi-linear parabolic equation (|6.4p satisfies 

\y{T)\^.^^^<Ce-^°^'^\yoy^^^ 
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Theorem 6.5. Let Q = (0,1), c satisfy and G satisfy ([631) - ([61]). There 

exists Co such that, for T > and M > 0, there exists positive constants C, ho such 
that for < h < ho and for all initial data yo € MP^ satisfying 12/01^2(^2) — there 
exists a semi-discrete control function v such that the solution to the semi-linear 
parabolic equation 

{dt-DcD)y + G{y) = l^v, ye(0,T) = 0, y(0) = yo (6.7) 

satisfies 

|y(r)L.(o) <Ce-'^°/"l2/oL2(n) 

where C = C{T,M). 

Observe that the constants are uniform with respect to discretization parameter 



A Proofs of Lemma 13.171 and intermediate results 
in Section [4] 

A.l Proof of Lemma [3.171 

We have 

(Cd£'u)„+3 - (cdZJu),,^! = hfn+l. 

As Du — pDv + Dpv we obtain 

''n+l(p„+|(c£'w)„+3 -p„+i (€£»«)„+ 1 + (i:'/9)„+3(cu)„+3 - {Dp),^^l{cv)„+lj 

=^h{rfU+,. (A.l) 
We write 

rn+lPn+l +rn+lPn+2 1 + (((T+)2p)r)„+i 

rn+iP„+3 = = := All, 

r„+i(crfDp)„^| = {rT+p)n+i{cdrDp)^^i 

= irT^P)n+i {{cdrdp)„+3 + {cdrDp)^^3 - {cdrdp)^^^^ 

= K21 (^{cdrdp)^^3 + , 

where K21 = (rr+p)„+i and K22 = icdrDp)„_^_i - {cdrdp)^^^. 
Similarly, 

Tn+lPn+l+rn+lPn 1 + (((T~)2p)r)„+i 

r„+ip„+i = = := iisi, 

rn+i{cdDp),^^i = {rT' p)n+i{cdrDp),^^i 

= {Tr^p)n+i [{cdrdp),,+ i + {cdrDp)^^i - (cdrdp)^^!^ 
:= K^i ((cdrap)„+ 1 + ^5:42'^ 
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where ^4:41 = (rr p)„+i and {cdrDp)^^j^i - (crfr9p)„+i . 
Additionally, 

(^)„+i = Vn+l + ^"^^ = Vn+i + 0{h){Dv)^+l, 

(^)„+| = Vn+l + ^"'^^ ^ ^"^^ = Vn+l + 0{h){Dv)n+l- 

From (jA.l[) we thus write 

i^ii(crfDt;)„+3 - K3i{cdDv)^^^i 

+ K2l[icrdp)n+s+K22){yn+l+Oih){Dv)„+s^ 

- X4i((crap)„+i + X42) + 0{h){Dvl^^^) 

= h{rf)n+i. 

Then 

i^ii((cd£'w)„+| - (cdDw)„+i) + (i^ii - K3i){cdDv)n+i 

+ K2l[*rCrdp]aV„+i + {K21 - K 41) {crdp)n+i Vn+l + {K21K22 - KiiKi2)Vn+l 

+ K2i[{crdp)n+^^+K22)0{h){Dv),,+ .+K,i({crdp)n+^^+K,2)0{h){^^^^ 
= h{rf)n+i- 
Moreover, as rdp — —Xs(j)dtp — sO\{l) we have 

w)„+3 - {cdDv)„^i^ + {Kii - K3i){cdDv)n+^ 

= K2lXs[-kC(l)d'4j]aVn+l - KVn+1 + (K2lOx{sh) + i^2li^22C(/l)) ({CdDv)n+^ " {cdDv)n+^^ 

+ (K2iOx{sh) + K2iK220{h) + K^iOxish) + K^iK^20{h)) {cdDv)^^i 
+ Hrf)n+i, 
where 

K = iK2i - Kii)icrdp)n+i + K21K22 - KiiK^2 
= (X21 - Kii).sOxil) + K21K22 - KiiKi2- 

From that, we can write 

L[{cdDv)n+-^ - (crfL>w)„+l) 
= K2l\s[^C<i)dlp]aVn+l - KVn+1 + H{cdDv)^^i + /i(r/)„+i, 

where 

L^Kii- K2iOx{sh) - K2iK220{h), 

K = {K21 - K4i)sOx{l) + K21K22 - K41K42, 

H = K2iOx{sh) + K2iK220{h) + KnOxish) + KiiKi20{h) - Kn + K31. 
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As L = 1 + Ox^si{sh) ^ (see below) then we read 

= (L-^K2iXs[kcc^di^]a - L-'K)vn+i + L'^ H (cdDv) 1 + L-^h{rf)„+i. 



We set 

Ji = L-^K2iXs[*c(l)d^]a - L-^K, 

J2 = L ^iJ, — L ^. 

We thus have 

(cd-Dt^)n+3 - {cdDv)^,^i = Jiw„+i + J2(cdL'w)„+i + J3h{rf)n+i (A.2) 

By using Proposition 13.101 we find 

l + (((T+)^p)r)„+i 
^11 = ^ = 1 + OxM^h), 

l + (((T-)V)r)„+i 
-fVsi = ^ = 1 + OxM'^h), 

K21 = {rT+p)„+i = 1 + OxAsh), 
K41 = {rT^p)n+i = 1 + Ox^sh), 
K22 ^ {cdrDp)^^3 - {cdr^p),^_^_3 = sOx,si{shf, 

K42 = {cdrDp)^,^i - (cdr9p)„+i = sC'a,j?(s/i)^ 

From that we estimate 

K = (K21 - Kii)sOx{l) + K21K22 - KiiKi2 = sOxAsh), 
H ^ K2iOx{sh) + K2iK220{h) + KnOxish) + KiiKi20{h) - K,i + i^ai 
= OxA^h), 

L = Kn - K2iOxish) + K2iK220{h) = 1 + OxAsh). 
For sh sufficiently small we have L^^ = 1 + OxA^^) ^^^^ then we obtain 
Ji = L^^ K2i\s[kc(t>dilj]a ~ L^^K 

\ + OxAsh))\s[^c(t)d^^]a + sOxAsh), 



J2=L-'H = OxAsh), 
J3 = = 1 + OxA^'^h). 

By using Proposition 13.141 Lemma 13.81 and Lemma 13.61 yield 
dtKn = 9t((((T+)V)0n+i) = Te{t)OxAsh), 
dtK^i = at((((T-)2p)r)„+i) = T9{t)OxAsh), 

dtK2i - at(rr+p)„+i = T9it)OxAsh), 
dtK^i = dt{rT-p)^+i = Te{t)OxAsh), 

dtK22 - dt({cdrDp)^+. - {cdrdp)^+.) = sTe{t)OxAsh)\ 
dtK^2 = dt({cdrDp)^+. - icdrdp)„+A = sT9{t)OxAsh)^ 
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which give 

dtL-^ = -j^ = (l + 0A,i?(s/i)) {dtKn + dtK2iOx{sh) + K2i{dts)Ox{h) 

+dtK2lK220{h) + dtK22K2lO{h)) 

where sh sufficiently small and 

dtH = dtK2iOxish) + K2i{dts)Ox{h) + dtK2iK220{h) + K2idtK220{h) 
+ dtK^iOxish) + Kiiidts)Oxih) + dtKiiKi20ih) + KiidtKi20{h) - dtK^i + dt 

It follows that we have 

dtJ2 = Te{t)OxAsh), dfh = Te{t)OxAsh). 

Furthermore, we can write (IA.2[) in the simple form 

= \s[^C(j)dil\aVn+l + XsOx,tx{sh)Vn+l 

+ OxAsh){cdDv)^+i + (l + OxAsh.)) /i(r/)„+i, 
which yields the conclusion. 

A. 2 Proof of Lemma 14.41 

By using Lemma 13.21 in each domain f2oi, f2o2 , we have 



111 = 2 / cr^p DpD{cdDv)Dv 

= 2 / cr'^p~Dj)D{cdDv)D^ + 2 / cr^pTJ^D{cdDv)D^ 

= 2 / cr^p~D^CdD{Dv)D^ +2 / cr^pU^{Dcd)(D^f 
■^Qoi -^Qoi 

+ 2 / cr'^p~D^CdD{Dv)D^ +2 / cr^pTJ^{Dcd)(D^f 
= J2 [ cr^P^CdD{Dvf + 2^2 [ cr^pTyp{Dcd)(D^ 



We then apply a discrete integration by parts fProposition I3.5|) in each domain 
f^oii ^02 with dfloi = {0, a} and dflo2 ~ {f^jl} for the first two terms and we 
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obtain 
hi - 



2 2 

J2 [ D{ccy~pWp){Dvf + j cr^Wp{Dcd)(D^f 

{c-cypWp){l){Dv)l^^^^. - [ [ccypWp){a+){Dv)l^, 

^ Jo ^ 

{ccypDp)ia')iDv)l [ {ccypD^)iO){Dv)l 

Jo ^ 

-Y. I D{ccyf,Wp){Dvf ^2Y, [ cr^^D^{Dcd)(D^f+Yn. 



where 



Jo ^ Jo ^ 

y}?^[ {ccypWp){a-){Dv)i^^- f {ccypWp)mDv)i 

Jo ^ Jo ^ 

Lemma A.l. (see Lemma B.3 in [BHLlOa]) Provided sh < ^ we have 

c,r^P^W^{Dcd^) = sA0,;O(l) + sO xAi^hf) ^ 

r^P^'Di^ = ndp, + sOxAishf) = + sOa,.«((s/i)'), 

r^piDpi = r^p^r^Dpi = (1 + Ox,fi{sh))riDpi. 



Moreover, by Lemnia [3.3l and ProDOsition l3 . 5l in each domain Oqi, ilo2 we obtain 

2 

[ s\(t)(D^f < [ sX(f{Dvf = j s\^{Dvf -^^BT,< [ s\^{Dvf 
Jn' Jn' Jn' 2 J^v 



1=1 



snice 



BTi = sX(l)(a)(Dv)l,i + sX(l)(a)(Dv)l > 

'''' 2 2 

BT2 - 5A(/)(l)(i?«)2 ^ .sXc^{a){Dv)l^s > 



and = + /i'^C'a(I) then we can write 



sX<j){DvY < sX(f){DvY + sXh'Ox{l){Dvy 



Similarly, we have 



Thus 



< s PxM'^hy UDv)^ < s pxAshY iDv) 



hi>- sX^{c^4>{i;'f)d{Dvf - + Fn, 
JQ' 
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where Xn = J^, i'ii{Dv)^ with i/n of the form sX(j>0{l) + sO\,fi(sh) and 

^11 - ^11 + ^11 + ^11 > 

Jo ^ 

T 

^1 + Ox,^{sh)){ccd)mrDpm{Dv)l, 











Jo ^ Jo 

A. 3 Proof of Lemma 14.51 

We set q — rpccf) . By using a discrete integrations by parts (Proposition 
and Lemma 13.21 in each domain fioi, ^02 we have 



2 r 

h2 = -2V / sqvD{cdDv) 

= 2^/ sqcd{Dvf + 2^ I sDqcdvDv 

sq{a^)v{a){cdDv)^_^_i + / sq{a+)v{a){cdDv)„,;i 



Jo 

2 I. 2 



2 2 

= sqcd{Dvy + 2^ sDqcdvDvdt + Yi- 



.l-^Qoi i = l-'Qoi 

since f |ano = and with 9rioi = {0, a}, 9r2o2 = {a, !}• 

Lemma A. 2. (see the proof as given in Lemma 4-4 of JEHLlOaf ) Let i = 1,2. 

Provided sh < ^ we have 

q, = r,p~ c0- = X^cM^'if + A0»O(1) + Ox.jiish)^, 
q, = X^{c<l)Mf)d + A0,O(1) + Oxj^{{shf + /i), 

Dgi = D(rip~ ) ccp- + (r^^) £'(€01') = Ox,si{^)- 

Note that the proof and the use of Lemma IA.2I are carried out in each domain 
f^oii ^^02 independently. 
It follows that 



12, 
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Then 



with 



T 

Yi2 = / s\^(j}{a)v{a)[c{il;'Y*CdDv]a 
Jo 

+ 5i2w(a)(ci:>w)„+3 +5i2w(a)(cL»w)„+i, 

where i5i2,^i2 are of form s{\(f){a)0{l) + Ox,ji{sh)^) and 

Jq'o Jq'o 

where 

j.i2 = sA0O(l)+sOA,j?(/i+(s/i)2). 

A. 4 Proof of Lemma 14.61 

We carry out a discrete integration by parts (Proposition I3.5P in each domain 
f^oi, ^^02 with dVloi — {0, a] and 91^02 = {a, 1} as follows 



/i3 = / rpD{cdDv)dtV + / rpD{cdDv)dtV 



D{rpdtv)cdDv — / D{rpdtv)cdDv 

+ [ {rp){a')dtv{a){cdDv\,+ i - / ir~p)iO)dtviO)icdDv) i 
Jo Jo 

+ {rp){^)dtv{l){cdDv),^^^j^3- {rp){a+)dtv{a){cdDv),^+3 



/ D{rpdtv)cdDv - / D{rpdtv)cdDv 

J QQ^ J Qm 



{rp){a^)dtv{a){cdDv)„^i - / {rp){a^)dtv{a){cdDv) 



^0 
2 „ 2 



= / D{rp)dtVCdDv I r^idtDv)cdDv+Yi3, 

Qi Q2 

by Lemma 13.21 and with 

Yi3= {rp){a^)dtv{a){cdDv)^+i - {rp){a+)dtv{a){cdDv)^_^_s., 



n+i 

'■ JO 



as v\qq,„ = 0. 



41 



By applying Proposition 13 . 13l in each domain floi, r2o2 we find 
On the one hand, we have 

2 2 

|Qi| < J2 [ s-^OxAsh){dtif + j sOx,si{sh){Dvf 
i=l JQoi i^i JQoi 

= f s~^OxAsh){dtvf + f sOx.si{sh){Dv)\ 

Jqo Jq'„ 

by {dtvf <{dtvf in each domain fJoi, ^^02 and X^Li !aoi (^«^)^ = /o,, 

On the other hand, by an integrations by parts w.r.t t we write as 



2 2 



We observe that for sh < ei(A) with ei(A) sufficiently small we have rp > 
by Proposition 13. 131 The sign of the term at t = T and t = are thus prescribed. 
Furthermore, Proposition 13.141 leads to dt{riPi) — T{sh)'^90x_fi{l), so that, for 
sh < R we obtain 

Q2>J2 [ T{shfeOxA^){DvY - Ca,^(1) I (Dv{T)f. 



Thus, 



with 



^13-/ {s{sh)+T{shf9)OxA^)iDvf+ s-^OxAsh){dtvf. 
Yis = / {rp){a~)dtv{a){cdDv),^^i - / {rp){a+)dtv{a){cdDv)„^3 
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A. 5 Proof of Lemma 14.71 

We set q = c^r^{DDp)Dp. Observing that Dv — Dv we get 



2 c^r'^{DDp)Dpi Dv + 2 c^r^{DDp)Dpi Dv 

JQni ' " ' JQa2 ' " ' 

qD{if+ [ qD{vf 

JQ02 

Dq{if - f Dq{vf 



Qoi Qo2 

Jo ^ Ja ^ 
Jo ^ Jo 

2 



'21 



= -E/ Dqv^ + Y.\l {Dq){Dvf+Y^- 
= -i:f mvf+jZ^I {Dq){Dvf +Y^lHYi\ 

by means of Proposition 13.51 Lemma |3.2[ Lemma 13.31 in each domain iloi, ^^02 
independently and where 



^21 — ^21 + ^21 — ^21 + ^21 + ^21 ' 

Jo ^ Jo 



1+ 



v\Q){Dq)i 



«'(l)(^9)„+™+| -^^ v\a){Dq)^^._ 
^ ^ v\a){Dq)^,^.-^ I v\a)[Dq)^^. 



2 Jo ' ^'"^^ 2 JO 



as v\dno = 0. 

We note that vi = !^{Dv)i, u„+„+3 = — ■|(Z?u)„_|_,,„^3 . On the one hand, 

by Proposition [3TTU] we have q = s^O\^si{l)rDp in each domain iloi, ^02- It follows 
that 



4'-'^ = r.'OA,i^(l)(ri?p)(l)(«)^+„+,+ r.^OA,r.(l)(ri?p)(0)(t))i 
Jo ^ Jo ^ 

^ OxAshf{r-D-p){l){Dv)l^^^, + r OxMshfirD^WKDvA. 
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On the other hand, by Proposition l3.15l Corollarv l3.9l we have q = —c^{s(t)X)^{ip')^- 
s^O\{l) + s^O\^fi{sh)^ in each domain fioi, ^02- We thus obtain 



A</.(a))=^(-(c2V'')(a-)(5)2^. +(cV'')(a+)(^)^„, 



— -'21 ^21 ' 

where 

Jo 

Lemma A. 3. (see Lemma B.8 in \BHL10a}j } Provided sh < R we have 

Note that the proof and the use of Lemma IA.3I are done in each domain r^oi i 
f2o2 separately. 



We then obtain 

.(2) 



h r ^, ^ h 





T 

Jo 

We thus write /21 

/2I > 3 / X'^s'^cj.^^i^'fivf - f fl2livf - f V2l{Dvf + Y2U 

jQ'n JQ'r. JQ'r. 



where 

M21 = (sA0)3O(l) + 5^0^ ,,(1) + s^OxA^hf, V2X = sOxA^hf. 

V ^(1,22) ^(2) 

J^21 - 121 + -'^^21 + ^21 + ^21 • 

A. 6 Proof of Lemma 14.81 

We set q = (?r{DDp)(j)" and by Lemma we have v = v + K^DDv/A in each 
domain fioij ^^02- It follows that 



Qo 



I22 = —2 / sqvv — 2 / sqvv 
' Qqi Qoi 



2 I sqv^ - I ^q{DDv)v 
2 [ sqv^ - [ ^q{DDv)v. 
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Applying a discrete integration by parts (Proposition 13.51) and Lemma 13.21 in 
each domain floi, ilo2 yield 



-2 ^ / sqv' + '^D{qv)Dv + F^^ 



4''=- r -^q{a-)v{a){Dv)^,^+ !!^q^a+)v{a){Dv)^^s, 







where 



as v\ono = 0. 

In each domain Qqi, $^02, we have (/) ~ Ca(1) and from Proposition 13.131 we 
have q = s^C'a,.,i(1) and Dq = s^C'a,.,i(1). We thus obtain 



Jo 

Lemma A. 4. (see Lemma B.9 and Lemma B.IO in WHLlOaf ) Provided sh < A 
we have 

h^DDq, = s{sh)Ox^A{l). 

Note that the proof and use of above Lemma [A. 41 are done in each domain ftoi, 
rio2 separately. 

Futhermore, we have (p" = A^(?/'')^'^ + A(/)C'(1) in each domain floi, fio2- It 
follows that 



sqi = SI 



in each domain iloij ^02- 
We thus write I22 as 



/22 = -2/ c2s3AV(V'')V+ / fl22V^+ I l^22{Dvf +Y22, 

'Q'o 
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where 



^22 — ^22 + ^22 ■ 

A. 7 Proof of Lemma 14.91 

By means of a discrete integration by parts (Proposition I3.5P in each domain 
f^oii f^02, we obtain 



/23 - V / cr{DDp)idtV 
= 2, / cr[DDp)dtv V 



_ = 1 JQ0^ 

^ £ {cr{DDp)mdMO)ii^ {cr{DDp))[a-)dtv{a)v^^._ 
^ ^ {cr{DDp)){a+)dtv{a)v^+^_-^j^ {cr{DDp)){l)dtv{l)v^^„,^^_ 

by Lemma 13.31 and where 

2 



2 ,2 r 

Q2 = V — / D{crDDp){Ddtv)v, 



1 



4'^ = (cr(^I?p))(a-)9tt;(a)f}„+i - | (cr(^I?p))(a+)a,t;(a)f}„+3 

as Stiilaao = 0. 

With an integrations by parts w.r.t t we have 

2 2 

By means of Proposition 13 . 1 3l and Lemma 13.71 in each domain floi, r2o2 we get 

nDDp, = s^Ox,si{l), 

and we further have 

Lemma A. 5. (see Lemma A.l in \BL12f ) 



dt{cnDDp,) = Ts'60x,si{y)- 
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Note that the proof and use of Lemma IA.5I are done in each domain fioi , ^02 
separately. 

It foUows that 



2 „ 2 



as \v\^ in each domain r^oi, ^^02- 



Moreover, we observe that /oo ^>-M^) ^ /oo ^a,^(1)^'^- Then, 

Ql = / rs20OA,^(l)^'' + / +OA,i5(l) ^u_, ). (A.3) 



We have 



By an integration by parts w.r.t t and Lemma 13.21 in each domain Hqi, rio2 we 



find 



02 = -Et/ dt{D{crDDp)i)Dv + Y.^ I D{crDDp)D{vf\lzS. 

* JQ0^ fcl ^ "'iloi 



By means of Lemma 13.21 and a discrete intergration by parts in space (Propo- 
sition 133]) in each domain JIqi, ^02 we see that 



= dtiDD{crDDp)y -Y^^ [ D{crDDp){dtv)Dv 



4=1 



2 '7 2 

^ v^{a)dtiD{crDDp)),,+ .+^ v^{a)dt{D{cTDDp))^^^ 







^ ^ / 9*(5i5(crDDp))t;2 -j^^f D{crDDp){dti)Dv + Y^^ 



as w|ano = 0. 

Lemma A. 6. (Lemma A. 2 in WL12f ) Provided sh < M. we have 

h^DD{c,r,{DDp,)) = s(s/i)C'a^^(1), 
h^dt{DD{c,nDDp,)) = Ts^eOx,f,{l), 
hdt{D{c,nDDp,)) = Ts'eOxA^), 
D{c,nDDpi) = s^Ox^l). 

Note that aU above terms are done in each domain iloi, ^^02 separately. 
We thus obtain ^ 

= [ sTeOx,A{shy{a). 
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Applying the Young's inequality and using that \dtv\^ <\dtv\^ in each domain 
f^oii ^^02, we have 

JQ'o JQo JQ'o 

(A.4) 



By using ProDOsition l3.5[ Lemma [A.Gl in each domain iloi, ^^02 separately yield 



,^1 » Jno. 



+^v\a){D{crDDp))^^^\lzl - ^v\a){D{crDDp)) ^^^^\lz^ 

sOxs.{sh){vf\t=T+ I sOxAsh){v)\=^ + OxAshfv''{a)\\zl 

sOxMsh){vf\t=T + I sOxA^h){v)\=^ + Yi\ (A.5) 

as v\dno — where 

Y,f^O,A^hrv'ia)\lzl 
Collecting (TOI) . (|X4|) and (|X5|) we obtain 

I23 > I s''{OxA^K + <^=^A^H J - ^23 + 1^23, 



where ^23 and I23 are as given in the statement of Lemma 14.91 

A. 8 Proof of Lemma [4.101 

By means of a discrete integration by parts (Proposition [33]) in each domain 
i^oi, ^^02 separately, we get 



'31 



= -2t / {dt9)ipcrDpvDv -2t / {dt9)ipcrDpvDv 

= -2t {dt9) ifcrTJ^v Dv - 2t / (dtO) ipcrU^v Dv + Y. 



with 



h r ... ^ ..^ . h '■^ 



as v\dno = 0. 



We have ipcrDpv—ipcrDp v + ^D{ipcrDp)Dv in each domain iloi, ^^02- It 
follows that 



48 



'31 = -rJ2 [ {dtO) {crDp^) D{vf -j^r^ I {dt9)D{crDp^){Dvf +Yjp 
r ^ / idte){D{crD-p^)y - E / idt9)D{crD^^)iDvr + Y^P 



Jo Jo 



,^iJQo. ^ JQ<H 

by using a discrete integration by parts in each domain fioi, f^02 separately and 



as w|af2 = 0. 

By using the Lipschitz continuity and Proposition 13.131 we get 



The proof is done in each domain iloi, ^^02 separately. Note that max9t^? = T6^ 
It thus follows that 

jQ'r. JQ'n 



where 



^31 - ^31 + J^31 ' 

Jo 
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A. 9 Proof of Lemma [4.131 

We see that 



^(1) ^(1,1) 

^11 + ^21 



^ (1 + 0,Msh)){cc,){l)irDpUDv)l^^^, 

T 







(1 + OxMsh)) (ccdmirDpUDv) 



2 



Jo ^ Jo ^ 



Moreover, by (|4.ip we have + i^2'"i — ^o^' sufBciently smah. 

We next focus our attention on the trace term at 'a' on Y-l^^'^'^ + Fj'i'^^'' as 
follows 

A. 10 Proof of Lemma 14.141 



f~\2 /Wn+l + t;„+2x2 / h ,2 

.^^)„+3 = ( ^ ) ^ {Vn+l + -{Dv)^+3) 



/j2 

+ —{Dv)l^3 + /iv„+i (£>«)„+ 3 

(A.6) 



Similarly, we have 



iv)l+l = + 4(erf)2 - (cdj:'«)„+i. (A.7) 
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We thus write 1^2^i follows: 

V(1'21) 

Jo 
Jo 

n+2 "'+2 

+ r(sA</.(a))3((c2v'^)(a+)— ^(c,Z)^)„+| - (cV'^)(a-)— ^(c,D^;)„+i)t;(a). 
io ^ 2(Cd)„+| = 2(cdj„+i 2/ 

(A.8) 

Moreover, the term ' is given by 

We estimate as 

c{a-)ci{a) ((cd)„+i + 0{K)) ((cd)„+i + (cd)„+|) 



((qUi+0(/i))(2(c,Ux+0W) 



Similarly, 



c(a+)c,(a) ^^^ 



(2 1) 

We thus obtain F^^^;^ ' 



n?''' = y^sA<^(a)(-^'(a-)(l + /iO(l))(cdi?<^i +V''(a+)(l + /iO(l))(cdZ?<+3) 
= rsA<^(a)[V'*(cdi?i;)2]„ 

JO 

+ / 8M{a)^\a^)0{h){cdDv)l^,_+ I 8\<^{a)i,\a-)0{h){^^^^ (A.9) 
Jo ^ Jo ^ 
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Combining (IA.8I) with (IA.9|) we obtain 



^11 "T ^21 

sA0(a)[^'*(cdi^^;)2]a+ / s3A3(/.3(a)[(^')'c'*]a«'+i 
Jo 

T i-T 



sA(/.(a)Vy(a+)0(/i)(c<ii?i;)2^3 + / s\<t>{a)^' {a~)0{h){cdDv)l^, 
^ Jo ^ 



where 



T /.T 

21 I / ^3\3i3/„\iJ2/j,/\3,] „,2 



Jo Jo 



n+l- 



and /zi can be written as 

Ml = / sOA,^(s/i)(cdi?^')' + 3 + / sOx,ii{sh){cdDv)l^^ 
Jo ^ Jo ^ 

Ox^Si{sh){cdDv) Vn+i + / s^Ox^si{sh){cdDv)^^iVn+i- 

"Jo 

We can write 

[{4>') * {cdDv)\ 

= [{'^')*]a{cdDv)l^i + [*{cdDv)]lij'{a+) + 2mcdDv)]aij'{a+){cdDv)„^ 
Indeed, we have 



[i^')*icdDv)X = {cdDv)l^si;'{a+) - {cdDv)l^,,p' (a-), 

and 



[{^')*UcdDv)l^^ + [^CdDv)]li^'{a+) + 2[-^icdDv)]a^/{a+){cdDv)^+i 
= {cdDv)l^,lb'ia+) - {cdDv)l^,J/[a-) 

+ {cdDv)l^,^^'{a+) + {cdDv)l^,^4>'{a+) - 2(crfi?i;)„+| (cdi?«)„+ i ^'(a+) 
+ 2{cdDv)^^3{cdDv)^+i^'{a+) - 2{cdDv)l^,^i^' {a+) 
= {cdDv)l^,J/{a+) - {cdDv)l^,J/{a-). 
Moreover, by using Lemma 13.171 we obtain 
[{^')*{cdDv)X 

= M)*UcdDv)l^^ + [*{cdDv)]ltiy{a+) + 2mcdDv)]ai^'{a+){cdDv)^^i 
= [{^P')*]aicdDv)l^, + (^\h^[^c^^')]lvl+^ +rl + 2Xsro[*c^^']aVn+i)^p'{a+) 
+2(^sX[*c4>^']aVn+i + ro^ip' {a'^){cdDv)„_^_i, 
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which gives 



II = i s\(j){a)[%l)'-k]a{cdDv) 
Jo 

+ I 2s2A2^(a)[*c0V']a/(a+)vjv+i(cd£>?;)„+i 
Jo 



Jo Jo 



+ 2 / sX(l}{a)il)'{a+)ro{cdDv)„ 
Jo 



>n+h. 

Moreover, we have: 

[^c<^V']a = [c^i^'*]l + 2[c(/.^'^]„/iOa(1) + h^Ox{l) = ^^a)[ci;'*]l + hOx{l). 
We thus write /i as 



fi = [ sX(l){a)[ilj'-k]a{cdDv)l 
Jo 

+ [ 2s^X''(l>''ia)[c^'*]a^'{a+)vn+i{cdDv)^^i 
Jo 

+ £ s3A3</.3(a)([cV'^]>'(a+) + [c'm'*]a)vl+, 

[ sX(l>{a)i^'{a+)rl+2 ( X^ (f" {a)[cip' *]aip' {a+)rQVn+i 

Jo Jo 



+ 



I-T 

+ 



2 / sA<^(a)V'(a+)ro(cdDt;)„+i + / s^Ox^^ishy^^^ 
Jo Jo 

+ / sOx,A{sh)Vn+l{CdDv)^j^i + i sOx,R{sh)rQVn+l 

Jo Jo 

= [ sX<t>{a)[i^'*UcdDv)l., 
Jo ^ 

+ [ 2s2A2</,2(a)[cV>'*]„^'(a+)v„+i(cd£'u)„+i 

^0 

X'<P'{a)([ci;'*]li;'{a+) + [c2(V')'*]a)t;^+i + Mr 
where can be written as 

lir = f sOxiiyo + I s'OxiiyoVn+i + / sOx{l)ro{cdDv)^+x 
Jo Jo Jo 

+ I s'^Ox,A{sh)vlj^-^ + / sOx,^{sh)vn+\{cdDv)^j^i + / sOx,A{sh)rov. 
Jo Jo ^0 



T 

+ / S^' 
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We have thus achieved 

/i = / sA(/>(a) (Au(t, a), a)) + /Xr, 

with u(t,a) — (^{cdDv)^^i, sX(l){a)vn+iy and the symmetric matrix A defined in 
Lemma 12.11 

From the choice made for the weight function /3 in Lemma l2. II we find that: 

/i > Cao / sX(j){a){cdDv)l^i + Cao [ (j)^ {a)vl+;^ + fir, 
Jo ^ Jo 

with ofQ > 0. 

A. 11 Proof of Lemma 14.151 

By using Lemma 13.171 we have 



^4 .0 

T 



T 



rp{a )dtv{a){cdDv)„^i, 







where Ji , J2 and J3 are given as in Lemma 13.171 

Since J2 — 0\^fi{sh) and rp = 1 + Ox,si{sh) we can write 

Yi3 = / Ox,si{sh)dtv{a){cdDv).^^i - / r/5{a+)Jiv{a)dtv{a) 
Jo Jo 

T 



rp{a^)J:idtv{a)h{rf)n+ 







Futhermore, as / = /i — dt{pv) we thus find 

Yia = / 0\^si{sh)dtv{a){cdDv)^,i - / rp{a^)Jiv{a)dtv{a) 



T 

rp{a+)J3dtv{a)h{rfi - rdt{pv)) 
With an integration by parts w.r.t t for the second term above we obtain 
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n+1 



^3 = / Ox,A{sh)dtv{a){cdDv)^^i+^f dt{rp{a+),hy{a) 
Jo ^ ^ Jo 

- ]^r-p{a+)J,v\a)\\zl - j r-p{a+)hdtv{a)h{rh)^+^ 

T 

rp{a'^)J3df v{a)hrn+i {pdtv + dtpv) 

OxAsh)dtv{a){cdDv)^+i +]- [ dt{rp{a+).hy{a) 

^ ^ Jo 

+ sOx,si{l)v^{a)\\z^ + ( Ox.n{l)dtv{a)h{rh)^+^ 

Jq 

+ [ {l + Ox.si{sh))h{dtv{a)f + ]- [ r^{a+)J3h{rdtp)n+idt{v^(a)), 
Jo ^ Jo 

where rp, J3 are of the form 1 + 0\,si{sh) and Ji of the form sO\^fi{l). 
We apply an integration by parts in time for the last term 

Yi3 = / OxAsh)dtv{a){cdDv)^+i+]-( dt{rf>{a+)J^y{a) 
Jo ^ ^ Jo 

+ sOx.si{iy{a)\\zl + I Ox.^{l)dtv{a)h{rf^U+^ 

Jo 

+ / {l + Ox,si{sh))h{dtv{a)f -]- I dt{rp{a+)h{rdtp)n+i)hv\a) 
Jo ^ Jo 

+ \rp{a^)Jj,{rdtp)n^ihv^{a, )\\% . 
Moreover, we have 

dts = s(2i - T)e = sTeO{l), 

dtp = -^{x){dts)p = -^{x)s{2t - T)ep, 

rdtp = -ip{x)s{2t - T)e (A.IO) 

dtirdtp) - sT''e^O{l), 

by using (E^- (lO) . 

Now we estimate the terms dt{rp{a'^)Ji) and dt{rp{a'^)J3{rdtp)n+i)- By re- 
calling dfJi = sT90\_si{sh), dtJz ~ T90\,si{sh) as well as using Proposition l3.14l 



and (|A.10|) we obtain 



dt{rP{a+)Ji) = dt{rp{a+))Ji+rf>{a+)dtJi 



and 



dt{rp{a+)Jz{rdtp) n+l I 
= dt{rp{a^)),h{r dtp) n+l + rp{a+)dtJ3{rdtp)n+i + rp(a+) Ja^t ((r9tp)„+i] 
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Thus Yi3 can be written 



Yi3 = / Ox,ii{sh)dtv{a){cdDv),^+i + / sTeOxM''^h)v'{a) 
Jo Jo 

,.2/„\\t=T 



+ sOx,si{lK{a)\\z'o + / Ox,si{l)dtv{a)h{rh)^+^ 

Jo 

+ / {l + OxA^h))h{dtv{a)f + / T^6^0x,s,{shy{a) 
Jo Jo 

+ \{l + OxA-<m)h[ - v{a)s{t){2t - T)9{ty{a, .)) . 

We observe that for < s/i < 63 (A) with 63 (A) sufficiently small we have 

= 1 + OxAsh) > 0. 

Additionally, f{x) < then the last term of Y13 are non-negative. From that, 
we estimate Y13 as follows 



Yi3> [ Cx,sih{dtv{a)f + I {sTeOxAsh)+T^e^OxAsh)yia) 
Jo Jo 

+ sOxA^yi^)\l=I+ [ OxAsh)dtv{a){cdDv\,+ i + [ OxA^)dtv{a)h{rh)n+i. 

Jo Jo 



A. 12 Proof of Lemma [4.161 

On the one hands, as / = /i — dt{pv) we write 

[rDn+i = (rfi)n+i - {rdt{pv))n+i 

= {rf i)n+i - [{rp)dtv + {rdtp)v) 

V / n+l 

= {rf l)n+l - {dtv)n+l - sT0Ox{l)Vn+l. 

We thus obtain 

|(r-/)„+i|' < C((r/i)2+i + {dtv)l+,+s^T^e^Ox{l)vl+,). (A.ll) 
On the other hands, 

[Pl*P2] = bl*]a(p2)„+i + Pl{a+)[*P2]a, (A.12) 

and we recall 

{cdDv)„^3 - (crfL)u)„_|_l = [kCdDv]a = Xs[-kC(l)1p']aVn+l + Tq, 

where tq is given in Lemma 13.171 as 

ro = sOx,ii{sh)vn+i + Ox,si{sh){cdDv)„_^_l + hOx,ii{l){rf)n+i, 
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We then have 

{cdDv)l^3 = {cdDv)l^i + [*CdDv]l + 2[*CdDv]a{cdDv)„+i 

= {cdDv)l^^i + X'^s'^[kc(l)i>']lvl_^_^ + 2As[*c(/)V']aroD„+i 
+2Xs[*c(l)^p']aVn+i{cdDv)^+i + 2ro (cd-Dw)„+ 1 . (A.13) 

and we compute 

+sOx^si{sh)'^{cdDv),^^i Vn+l + sOx^ii{sh)h{rf)n+iVn+i 
+Ox^Si{sh){cdDv)^^ih{rf)n+i. 

By applying Cauchy-Schwartz inequahty we have 

{cdDv)l^, < 0{l){cdDv)l^, + s^Ox{^)vl+, + 0{l)rl (A.14) 

rl < s^OxA^hfvl^, + OxA''^hf{cdDv)l^^ + h'OxA^){rf)l+i, (A.15) 



srovn+i < {s^OxAsh) + sOxA^)H+i + OxAsh){cdDv)l^,+hOxAsh)irf)l+u 

(A.16) 

s^rovn+i < {s''OxAsh)+es^OxA^))^l+l+^^OK^ish){cdDv)l^^+C,hOxMsh){rf)^^^^ 

(A.17) 

sT0roVn+i < {s^TeOx.Ash) + sT^e''OxAsh))vl+^ 

+ sOxA^h)(cdDv)l^^^+hOxAsh){Tf)l+^, (A.18) 

sro(ci?«)„+i <s^OxAsh)vl+^ + {sOxAsh)+^sOxA^)){cdDv)l_^^+C,OxAsh)h{rf)l+^, 

(A.19) 

{dtv{a))ro < OxA^Mdtv{a)f+sOxAsh){cdDv)l^^+s^OxAsh)vl+^+OxA^)Hrf)l^^^ 

(A.20) 

We estimate following terms 

The first term, by using (|A.14p we have 

T 

sOxAsh)^icdDv)l^^ + sOxAshf{cdDv)l^, 



-"^11 





< / sOxAshf{cdDv)l^^+[ s^OxAsh)^yl+i+ f sOxA^hfrl 
Jo ^ Jo Jo 
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Moreover, by using (jA.151) we obtain 



T 



Then, by using (|XTT]) we estimate f/^'^' 



- 11 



Jo ^ Jo 



For the next term, using (IA.12[) and Lemma 13.171 we obtain 



T 



Yi2 = I sOx{lHa)[ci^'^*{cdDv)]a 



sOxA^Ha){cdDv)„+i + sOxA^Ho-){{cdDv)„^i + sOx{l)via) + rg 
sOxM'^Ha)icdDv),,+ i + s^OxMC^yia) + sOxM'^M'^yo- 



Using (fOel) yields 

sOx^Si{^)v{a)ro 

< (^s'OxMsh) + sOxA^))^l+i+ OxAsh){cdDv)l^^+j^ hOxAsh){rf)l+^. 

By using (jA.ll|) we obtain 

11^121 < r {s^OxA^) + sT^S^OxAsh)^yn+i+ r Ox.Ash)h{dtv)l^, 
Jo Jo 

+ f OxA^){cdDv)l^^ + f OxAsh)h{rh)l^,. 
Jo ^ Jo 

Moreover, we have 

W„+i = Vn+l - -{cdDv),^^l = Vn+l + ©(/l) (cd l , 

= Vn+l + 0{h){cdDv)^^._. (A.21) 
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By using (|A.21j) . (|A.14j) we obtain 

{v)l+^^ + (.v)l+^^ < 0{l)vl^, + 0[h?){cDv)l^,^+0{h^){cDv)l^,^ 

< (0(1) + 0{shf)vl^, + 0{h')icDv)l^^ + 0{h'y, 

Thus, we have the following estimate 



(1,22) 



< 



Futhermore, using (|A.15|) we have 



'0 "'0 ^ "'0 



By using (|A.11[) we get 



Y, 



(1,22) 



< 







For the term 1^22^'' "we have 



Y, 



(1) 



^ ;j2 

T 

T 

^fi{sh)^{cdDv)^^i + s'^Ox^ii{sh)^v{a) + sO\^ji{sh)^ro^v{a) 



Using (IA.16|) we achieve 

i-T 

sOxA^h)'^via)ro 

< 1^ s^OxAshfvl+i+ \^ OxAshf{cdDv)l^^+ 1^ hOxAshf{rf)l^,. 
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Using (IA.11[) . we estimate 1^22^'' ^s: 



^22 



(1) 



' 22 



2 

n+l 



And, using (jA.2ip and Lemma [3. 171 we obtain 







sOxAsh){dtv{a))v{a) + f Ox^sh)^ {dtv{a)){cdDv) 
"'0 







T i-T 

2/ 



sOx,si{.sh){dtv{a))v[a) + \ OxAshr{dMa)){cdDv),,+ i 
Jo 



In addition, with s, A enough large, sh enough smah and with applying Young's 
inequality and (|A.20p yield 

sOxMsh){dtv{a))v{a) < [ OxAsh)h{dtv{a)f + [ s'OxA^yia)- 
Jo Jo 

OxAshfidtv{ajKcdDv\,+ i< [ OxAsh)h{dtv{a)f+( sOxAshf{cdDv)l,^. 



OxAshr{dtvia))ro < / OxAshyh{dtv{a)y + / sOxAshricdDvy^, 
Jo Jo ^ 



T i-T 

3/'„a2 



+ / s'OxA^hY{v)i+, + / OxAshYh{rf)i^^. 
Jo Jo 



Using (|A.11[) . we estimate ^23^'' 



yd) < 
^23 — 



[s^OxA^) + sT^O^OxAshf) vl+, + j OxAsh)h{dtv{a)) 



^ 'Jo 

+ / sOxAshf{cdDv)l^, + / OxAshfKrh)l+i- 
Jo ^ Jo 
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Applying Lemma [3. 171 we have 



T 

sr6'C'A,.«(s/i)w(a) (c<j 1 

sr0C>A..«(s/i)w(a) (sOA,.s(l)w(a) + 
By using (|A.18p we obtain 
^ sT0OA,^(s/i)v„+iro < ^ (s2r0OA^.«(s/^)2 + sT2020^^^(^;j)2^^: 



2 

n+1 



Jo ^ Jo 

We have 

sTeOxAs'T'){cdDv)^+iVn+i < [ sT^e^Ox,si{sh)vl+,+ [ sOx.si{sh){cdDv)l^r. 

Jo Jo ^ 

(1) 



With (lA.lip we thus estimate Yg^ ' as 



T 

(1) 



^31 — 

JO ^ 'JO 

Jo ^ Jo 

Next, by using (jA.14l) we estimate /i^^' as 

Jo ^ Jo ^ 

< f sOxAsh)icaDv)l^, + [ s^OxAsh)vl+i + [ sOxs^ishy^ 
Jo ^ Jo Jo 

By making use of (jA.lsp we have 

i-T 

sOxAshyl 



< 





^ s'OxA^hfvl+i + r sOxAshf{cdDv)l^r + C OxAshfKrf)l+i- 
Jo ^ Jo 



Using (|XTT|) we obtain 



(^s^Ox,dsh) + sT^e^OxAshf) ^l+i + / OxA-^hfKdtv)l+^ 



+ r sOxAsh){cdDv)l^^+ r OxAshfKrhfn+i- 
Jo ^ Jo 
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By making use Lemma 13.171 we have 

Jo Jo 

s'^Ox,M{sh)vn+i{cdDv)„_^i + / s^Ox^M{sh)vl_^i+ / s'^Ox^si{sh)rov 
Jo Jo 

Applying Young's inequality and using (|A.17[) yield 
^ Jo Jo ^ 



T 

2 



s Ox.,iiish)roVn+i 







< r s^OxAshfyl+i+ r sOxMsh)HcdDv)l^r + C OxA^hfKrffn+i- 
Jo Jo ^ Jo 

Using (TOT]) we have 



Mi^' < 



(s^OxAsh) + sT^9^0xAshf)^l+i + / OxMshfHdtv)l+, 



T 



+ / sOx,MishKcdDv)i_^, + / OxMshrhirh)i+,. 
Jo ^ Jo 

We thus obtain 

m < ^ (s''OxAsh) + sT^e^OxAshf)<+i+ OxAshfHdtv)l+, 

+ r sOxAsh){cdDv)l^, + r OxAshfh{rh)l+,. 
Jo ^ Jo 

Now, we estimate some terms of fir- By using (jA.lsp - (|A.19p we have 



sOx{l)i 











< 



^ s'OxMshfvl^,+ r sOxA^hf{cdDv)l^^ + f OxAsh)h{rf)l^^. 
Jo ^ Jo 



s'OxiiyoVn+i < ^ (^s^'OxM'^h) + es^'OxA^))^ 



2 

n+1 



+ r sOxAsh){cdDv)l^,+C, r OxAsh)h{rf)l^^. 
Jo ^ Jo 
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Jo Jo 

i-T 



1 . 



sOx,si{sh)vn+i{cdDv)j^^i < / sOx,si{sh){cdDv) 1 + / sOxAsh)v^_^_i. 

Jo ^ Jo 

Using (IA.11|) we have: 



+ I 0,,xAsh)h{dtv)l^^ + / 0,,xA{sh)h{rh)l^^ 







+ ^ {sOxAsh)+esOxM^))icdDv)l^,_. 
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